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ABSTRACT 

We develop a convergence theory for semi-discrete approximations to nonlinear systems 
of conservation laws. We show, by a series of scalar counterexamples, that consistency with 
the conservation law alone does not guarantee convergence. Instead, we introduce a no- 
tion of consistency which takes into account both the conservation law and its augmenting 
entropy condition. In this context, we conclude that consistency and L°°-stability guar- 
antee for a “relevant” class of admissible entropy functions, that their entropy production 
rate belong to a compact subset of t). One can use now compensated compactness 

arguments in order to turn this conclusion into a convergence proof. The current state 
of the art for these arguments includes the scalar and a wide class of 2 x 2 systems of 
conservation laws. 

We study the general framework of the vanishing viscosity method as an effective way to 
meet our consistency and L°°-stability requirements. We show how this method is utilized 
to enforce consistency and L°°-stability for scalar conservation laws. In this context, we 
prove under the appropriate assumptions (L°°-bounds), the convergence of finite-difference 
approximations (e.g., the high-resolution TVD and UNO methods), finite-element approx- 
imations (e.g., the Streamline-Diffusion methods) and spectral and pseudospectral approx- 
imations (e.g., the Spectral Viscosity methods). 
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1. INTRODUCTION 


We consider 27r-periodic initial- value problems, which consist of the one-dimensional 
system of conservation laws 


(l.la) 


1)] + A[/(t i(z, i))] = 0, 


together with an augmenting entropy condition, which requires that for aU convex entropy 
pairs ( U(u),F{u ) = /“ < U'(w), f'(w)dw >) the following entropy inequality holds [7] 


(1.16) 


^[f/( U (x,!))l + ^[F(u(x,i))l<0. 


We want to solve this problem by a semi-discrete algorithm. To this end one associates 
with a large parameter N (or a small parameter Ax = , depending on the point of 

view) , i 

I. A finite (2N + l) -dimensional space, </#, spanned by 2^-periodic basis functions, 

and 

II. A possibly nonlinear, /-dependent, spatial discretization operator, P ! N , 


P f N : L°°[ 0,2 tt] -> $ N , 

such that Pjyio(x) is an appropriate $ ^-approximation of f(w(x)); here, denotes the 
(2 N + l)-space spanned by the 27r-periodic primitives of {^(x)}^.^. 

The exact solution of (l.la), (1.1b), u(x, t), is then approximated by a ^-element, ujv(x, t), 

N 

(1.2) u N (x,t) = Y, MtfPkix), 

k=-N 

which is determined by the following procedure: 

Starting with prescribed (possibly pre-processed) initial-data u^(x, 0) in </#, we let 
ux(x,t) evolve later in time according to the (2N + l)-dimensional approximate model 

(1.3) *)1 + §^ P n u n{x, f)] = 0. 

In what sense does (1.3) approximate (l.la)? Let us rewrite (1.3) in the form 

d d d 

(1.4) [u N (x, *)] + — [/(u^(z, *))] = — [f{u N [x, t )) - P f N u N {x , *)]. 

The expression inside the right brackets is the discrete local error, E { N u^{x,t), 

(1.5) E f N u N {x,t ) = f{u N [x,t)) - P{jU N (x,t), 
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which reflects the amount by which u^(x,t) fails to satisfy (1.1a). Its size may serve us as 
a measure for the order of accuracy of (1.3). 


Definition 1.1: ( Order of accuracy) The approximation (1.3) is accurate of order 
■s > 0 with the conservation law (1.1a), if there exist constants, Cj,j = 0, l,***,s — 1, 
such that for any r,0 < r < s, and for ah wj^(x) in fa, the following estimate holds 1 

Co ‘ II + ]C c y • 

ox ,=x 

Here Cj,j = 0, 1, •••, are constants which may depend on the Z/^-bounds of w#(x) and 
its first j-derivatives, but otherwise are independent of N. 

Remark: In a similar manner one can define accuracy of nonintegral orders s > 0. 

We note that (1.6) is a refinement of the usual definition of accuracy, in that here, the 
local error on the left and the highest derivative involved on the right are weighted by the 
L 2 rather than the usual L°°-norm. This then implies that s-order accurate approximations 
are, in particular, consistent with the conservation law ( 1 . 1 a), in the sense of 

Definition 1.2: ( Consistency with the conservation law) The approximation (1.3) is 
consistent with conservation law ( 1 . 1 a) , if for all Wjv(x) in 4>n which are uniformly bounded 
w.r.t. N, there exists a vanishing sequence, cn — » 0, such that the following estimate holds 

(1.7) \\E{fW N (x ) || < c N \\^~w N {x)\\ . 

Indeed, an s-order accurate approximation satisfies the consistency requirements, (1.7), 
with Cjf = £“•. 

Remark : In the generic case, the discrete local error is upper-bounded by 

(1.8) II^»WJ»(*)II < ^ L \\§ i f{^))\\ . 

with some universal constant Const. Hence (1.7) holds with c ^ = Const 0 • where 
Consto ~ || A = f'(w N (x )) ||£°°(*). 

The first of the two ingredients involved in discretizing the conservation law (l.la), is 
the choice of the finite (2iV + 1 ) -space < 5 6 ^. A very convenient choice in this respect is 

1 We use < •, • > and | • | to denote the usual Euclidean vector inner product and norm. Similar notations 
are used for vector functions, e.g., (-, •) = / 0 2 " < > dx and || • || 2 = (■, •) for spatial vector functions and 
(v) = fo So " < V > dxdt, || • || 2 = (•>•) for Lf oc {x,t) space-time integration. 


(1-6) \\EfflVff(x)\\ < jp 


2 





the space of N - trigonometric polynomials, ttjv, spanned by In this case, 

coincides with the space spanned by its 2^-periodic primitives, and they consist of spatially 
smooth approximants of the form 

(1.9) u N (x,t) = ]T u k (t)e' kx . 

k=-N 

The principal raison d’etre for the tt# space is contained in its isometry with the (2JV + 1)- 
space of Fourier coefficients {u k (t)}^ = _ N . Spectral and pseudospectral methods are the 
canonical examples for TT^-discretizations of the conservation law (1.1a). 

Example 1.3: Spectral Methods. Denote by S^w(x) the spectral-Fourier projection 
of w(x) into ft pj 

(1.10) Sn w(x) = V w(k)e ,kx , w(k) = — - f w(x)e~ tkz dx. 

k=-N 2?r Jo 

The approximation (1.3) with P^u^(-,t) = S N f(u N (-,t)), amounts to the spectral method 
for the trigonometric approximant, u N (x,t) = u k (t)e ikx , which reads 

(l-H) — Ufc(t) +ikfk{t) = 0, f k (t ) = — y o f{u N [x,t))e" ikx dx, |fc| < N. 

In Section 7 we prove that the spectral method (1.11) is accurate of any order s > 0, i.e., 
it is “infinitely-order” or spectrally accurate with the conservation law (1.1a). 

Example 1.4: Pseudospectral Methods . Denote by rpNi^ix) the t/’dospectral-Fourier 
projection of to(x), which interpolates iu(x) at the 2 N + 1 equidistant collocation points 
x v+ 0 = (v + 0) Ax, v = 0, 1 , • • • , 2 N, with fixed 0 < 0 < 1 and Ax = irf + i , 

N a t 2N 

(1.12) il) N w(x) = ™{k)e' kx , &(/:) = —^ty(x v+ «)e 

k=-N v =0 

The approximation (1.3) with P^u^(-,t) = ^/(u^ (•, i)), amounts to the pseudospectral 
method for the trigonometric approximant, u N (x,t) = ^ 2 k= ^ N u k {t)e ikx , which reads 

d A a_ 2 N 

(1.13) —u k (t) + ikf k (t) = 0, f k {t ) = — J2 f{ u N{x v + 6 ,t))e- ikx '’+\ |fc| < N. 

dt W i/=0 

In Section 7 we prove that the t/>dospectral method (1.13) is accurate of any order s > 
i.e., it is “infinitely-order” or spectrally accurate with the conservation law (1.1a). 

Finite-difference and finite-element methods are usually interpreted as evolution dis- 
cretizations of (l.la) in terms of piecewise-constant or piecewise-linear approximants. As 
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we shall now see, these methods can be equally well interpreted as ^//-discretizations of 
(1.1a), governing the evolution of trigonometric approximants instead of piecewise-constant 
or piecewise-linear ones. Let us turn to 

Conservative Methods . We deal with discrete approximations of (1.1a) which admit 
the conservative form [9] 

(1.14) ^u N (x v ,t) + ~ = °’ x„ = i/Ax, v - 0,l,---,2iV. 

Here h f , = h / (u//(xj,_ p+ i,t),*-*,u/y(x,, +p ,t)) is the Lipschitz continuous numerical flux 
which is consistent with the differential one 


(1.15) 


h f (u N (x,t),- • • ,u N (x,t)) = f{u N (x,t)). 


In Section 7, we prove that arbitrary conservative schemes of this form are at least 
first-order accurate approximations of the conservation law (1.1a). In order to interpret 
such schemes within the 7T//-framework (1.3) (1.9), let us introduce the ^//-polynomial, 
Hu{x, t ) = H k (t)e ,kx , which interpolates the numerical flux values at the grid points 


x 


v+ 


i = {v + |) Ax, u = 0, 1, • ■ • , 2 N, 


(1.16a) 


A-r 2JV 

H N {x = x v+ i_,t) = h f u+1 _, H k {t) = — ^I + i« ’ Xv+ K 


u=0 


and then define Pj/M//(x, t ) as the sliding average of this interpolant 
(1.166) P f N u N {x, t) = H H^(x) = ~ J2Z_ H N {i)di. 

By applying such averaging plus summation by parts to Hfj(x,t), we end up with the 
(possibly nonlinear) spatial operator, P^, which approximates f(u N (x,t)) from 7r N via 

N ■, 2N 

(1.17) P f N u N {x,t) = £ Pk(t)e ikx , p k {t) = T-rriZ[ h i+i. ~ h l- i} e ~ <kx ‘' ■ 

k =-N v=0 2 2 

We note that due to conservation, Po(t), is well-defined as po(t) = 0. The approximation 
(1.3) with Pjy defined in (1.17) amounts to 

. 1 2N N 

(1.18) — u k (t) + —J^[h f L - h[_ = 0, u N (x v ,t ) = u k {t)e ikz % 

at ^ i/=0 2 2 k=-N 

which by the inverse discrete Fourier transform is equivalent with (1.14). 

In this manner one can imbed any finite-difference or finite-element method as an evolution 
scheme in 7T//. 
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An instructive example for this process is provided by 


Example 1.5: Centered Finite-Difference Methods . Consider the standard centered 
finite-difference method 

(1.19) ~u N (x r ,t) + jlj[/(u w (x„ + i,f)) - /( Uri(x„-i,t))\ = 0. 

The discrete Fourier coefficients, pr(t) of P^u^(x t £) in (1.17), equal 

MO = - /(«»(**-!, 0)1«-"*" = 

( 1 . 20 ) 

= /(«„(*., 0)e-“". 

In Section 7 we prove that the finite-difference method (1.19) satisfies the expected second- 
order accuracy requirement (1.6). We note that the numerical flux associated with (1.19), 

(1.2l) = ^ /( ti ^( 2 'i / + 1 > 0)]> 

depends linearly on the gridvalues of /(u^(-, £)); consequently we end up with a corre- 
sponding discretization operator, P^, which operates linearly on /(«//(•, t)), i.e., 


(1.22a) 


Pn u n( •>*) = PD Ax f(u N (-,t)), 


where 


(1.226) 


FD Ax w(x) 


E *(*)«“*. 


k=-N 


w(k) = 


sinkAx 

kAx 


Ax 2N 

Z7r V - 0 


—ikx v 


The fact that the finite-difference operator, FD Axt turns out to be a smoothed version of 
the V'dospectral projection rp N in (1.12), 

1 fx+Ax 

(1.22c) FD Ax w[x) = — — / il>N{w(t))d£ 

l Ax J x—Ax 

is typical to all standard linear differencing methods, consult [19]. In other cases, the 
resulting operator Pfq in (1.17) may depend on the flux f in a more intricate way as shown 
by 


Example 1.6: Finite-Element Methods . Let (U*(u) , F t (u)) be any preferred entropy 
pair associated with (1.1a). Define the entropy variables 

(1.23) v = lf'(u), 

and note that the strict convexity of £/*(ti) enables to uniquely invert (1.23), u = u(v). 
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The finite-element approximation of (1.1a), based on piecewise-linear elements in the 
entropy variables v reads [18, Section 5] [21, Section 4], 

(1.24a) — ~tij^(x v _i,t) + — u/y(x,/,i) 4- — un^Xv+i, t)j + = ® > 

and the mass lumped version of this yields 

(1.246) 4- — == 

Here h* x is the numerical flux given by 
* /+ 2 

(1.25a) K + ,=£ ,/(«(»„,({; ))J(, 

" > — 2 

where 

(1.256) v„ + i(£) = [\ ~ (*«/,*)) + (^+ £)Ul{u N {x v+1 ,t)). 

We recall that in the case of standard linear finite-difference methods, the numerical flux 
depends solely on the gridvalues f(uN(x„,t)). In contrast, the current finite-element nu- 
merical flux in (1.25a), (1.25b), h* u+x , depend on all the intermediate values /(u(u(f))). 
This fact is reflected in the corresponding discretization operator, P^, which according to 
(1.16a), (1.16b) is given by 

(1.26a) P f N u N (->t) = FE Ax f(u N (-,t)) , 

with 

( 1 . 266 ) FE^w(x) = ± £ + J Z’ N (t)dt , B'„(x = x„ + , ) = /££ u.(t>„ + i(«)<ie . 

Remark : In both cases of the finite-difference and finite-element methods (1.22), (1.26), 
the spatial discretization, P^, operates linearly on the flux /, i.e., 

(1.27) P£ l/lWa = aiP h + a2 ph m 

The more modern shock-capturing techniques employ numerical fluxes, h f L , which depend 
on the data, Un(x, t) and /, in an essentially nonlinear manner [3]. These methods can be 
also interpreted as evolution schemes in 7 Tjv using the discretization recipe (1.16a), (1.16b); 
in such essentially nonlinear cases, however, the resulting discretization operator, Pjy, need 
not satisfy the linearity property (1.27). Our foregoing discussion equally applies to linear 
methods as well as the essentially nonlinear ones. 
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The examples considered so far, were formulated within the trigonometric 7Tjv-framework. 
We close this section with an example which shows how the finite-difference/element meth- 
ods (1.24) , can be equivalently formulated within the canonical piecewise-linear framework. 


Example 1.7: Finite-Difference/Element Methods revisited. Let be the space of 
piecewise-linear grid functions, u>w(x), 


(1.28 a) 


N 


w N (x) = £ tz; w (x*)A fc (x) , 

k——N 


spanned by the 27r-periodic ‘hat’ functions, tp k (x) = Ajt(x), which are centered at x k = kAx, 

(1.286) A fc (x) = -^-min([x] - x k . 1 ,x k+1 - [x]) , [x] = x[mod2^] . 

This family of basis functions is not orthogonal - the corresponding mass matrix M, 
Mij = (A,(x), Ay(x)), is given by the circulant matrix 


(1.29) 


M — AxC , C = Circulant ^ . 

\6 6 6 / 


In order to interpret (1.24) within the current piecewise-linear framework, we introduce 
the finite-element discretization operator, FE Ax , which maps L°°[— 7T, tt] into the space $j\r 
of piecewise-quadratic B-spline gridfunctions. 


(1.30a) 


N 


FE Ax w{x) = £ za t+ in fc+ i(x). 


k=-N 


Here, n fc+ i (x) stands for the 27r-periodic “bell-shaped” functions supported on [xj fc _ 1 ,x fc+2 ], 


(1.306) 


^jfc+iO C ) _ ■A-fc+i( a: )] J 


and w = {w k +i} k= -N is the vector of cell averages 


(1.30c) 


W k+ L = W 


^o s -LIT w[x)dx - 


Differentiating (1.30a) we obtain, in view of (1.30b) and (1.29), 


(1.31) 


N 


—FE Ax w{x) = £ — [w fc+ i-w fc _i]A fc (x) . 


k=-N 


Ax 1 


Hence approximation (1.3) with PlfU N (-,t) = FE Ax f[u N {-,t)) reads, at the collocation 
points x = x„, 


(1.32a) 


±u N (x l/ ,t) + -±-[f i] = 0; 
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since uj v(*,i) is piecewise-linear we have 

— 1 r x v+i fi/2 

(1.326) f v+ i = -^J x f{u N {x,t))dx = J (= _ 1/2 fi u v+kti)) d Z = K+k > 

and we recover the mass lumped version of the finite-element method (1.24b), with 17, (u) = 
|u 2 as the “preferred” entropy function, 

(1-33) — Ujv(a;„,f) + ~ = ® • 

Similarly, the finite-element method (1.24b) with no mass lumping corresponds to 

(1.34a) FE Ax w{x) = ™k+i n k+U x ) . 

k=~M 

where Aw = {w k+ i — u>k-~)kL-N i s a vector solution of 2 

(1.346) C ■ Aw = Aw , Aw = {u7 Jfc+ i — w k _i}^ = _ N . 

In this case, (1.31) is replaced by 

(1.35) ^-FE Ax w(x)= -^[C' 1 ■ Aw] k A k (x) , 

and approximation (1.3) with P^u^(*,t) = FE Ax f{uN{’,t)) recovers the finite-element 
method (1.24a) 

(1.36) — ^~tijy(£i/_i, f) + — t ijv(^j/>t) “I - ~^n{ x u+1: t)| + ’^^l^'u+k ^v_l] = ® ’ 


2 Since the vector of ones is an eigenvector of C, a conservative solution of (1.34b) exists. 



2. THE NECESSITY OF ENTROPY DISSIPATION 


The consistency condition (1.7) guarantees that as N tends to infinity, the approxima- 
tion (1.3) approaches the conservation law (l.la). The question whether the approximate 
solution, ujv(x,t), approaches the corresponding conservative solution, u(x, t), is the ques- 
tion of convergence . The following scalar counterexamples-one for each of the discretiza- 
tions methods mentioned above, show that the solutions of consistent approximations may 
fail to converge to the appropriate conservative solutions. 


Counterexample 2.1: [20]. The spectral-Fourier approximation (1.11) of the scalar 
equation (l.la) reads 


( 2 . 1 ) 


d d 

— [u w (x,i)] + ~[S N f(u N (x,t))\ = 0. 


Multiplying this by u^{x, t) and integrating over the 27r-period, we obtain that u^-being 
orthogonal to 4^[{I — Sn) f i u N[x,t))], satisfies 


( 2 . 2 ) 
i d r 2r 


2dtJo u Ni x ’ t ) dx = ~ j 0 u N{x,t)—[f{u N (x,t))]dx = ~J u/'(u)<£u|*zg x = 0. 

Thus, the total quadratic entropy, U(u) = |u 2 , is globally conserved in time 

1 t 2 * 2JV 2N 1 

(2.3) - / u 2 N (x,t)dx = X] U(u N (x u ,t)) Ax = Y, U(u^(x v ,0))Ax, U(u) = -u 2 , 

* J 0 ,,—n n * 


which in turn yields the existence of a weak L 2 (x)-limit u(x, t) = w lim^oo un(x, t ). Yet, 
u(x, t) cannot be the entropy solution of a genuinely nonlinear (GNL) Equation (l.la) 
where /"(•) ^ 0. Otherwise, <Sjv/(tijv(x,t)) and therefore f(u^(x, t)) should tend weakly 
to f(u(x,t)); consequently, since f(u) is GNL, u(x,t) = slimtf_ 00 u/v(x,t) which by (2.3) 
should satisfy | u 2 (x,t)dx = | f£ x u 2 (x,0)dx. But this is incompatible with the entropy 
condition (1.1b) if u(x, t) contains shock discontinuities. 


Counterexample 2.2: The pseudospectral approximation (1.13) of the scalar conser- 
vation law 


(2.4) 
reads 

(2.5) 


du d(e u ) 

1 — - = 0 , 

dt dx 


^MM)] + ^ K e“" W) ] = 0. 
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Multiply this by rJjjfe UN ^ x,t ^ and integrate over the 27r-period: since the trapezoidal rule 
is exact with integration of the 2iV-trigonometric polynomial obtained from the second 
brackets, we have 

(2 . 6 ) £ g = 0. 

Thus, the total exponential entropy, U{u) = e u , is globally conserved in time 

2 N 2 N 

(2.7) Y! U(u N (x u ,t))Ax = U[un[x u , 0))Ax, C/(u) = e u . 

u =0 i /=0 

Hence, if ux[x, t ) converges (even weakly) to a discontinuous weak solution u(x, t) of (l.la) , 
then rfj^e UN ^ x,i ^ tends (at least weakly) to e u ^ x,t K Consequently, (2.7) would imply the global 
entropy conservation of / 0 2 t e u ^ x,t ^dx in time, which rules out the possibility of u(x,t) being 
the unique entropy solution of (l.la). 


Counterexample 2.3: [8]. The centered-difference approximation (1.19) of the scalar 
conservation law (2.4), yields the nonlinear completely integrable system [5], [12] 

(2.8) ~u N {x v , t) + - e u w(*^-i.O] = o. 

at 2Ax 

Multiplying (2.8) by e UN ^ Xu,t ^ we obtain 

(2.9a) ^[e“ w(W) ] + ~ = °* 

where the consistent entropy flux, F v+ i = F(u N (x v ,t),u N [x l/+ i,t)), is given by 

(2.96) F + 1 = — gl u w(**'»0+«w(* , '+ i i t )l. 


In particular, this implies a global entropy conservation of the exponential entropy, 
U(u) = e°, 


( 2 . 10 ) 


2 N 2N 

) ] 17(ujv(x l> , t)) Ax = "y ' U(un{x„,0))Ax, 


u=0 


v=0 


U(u) = e u . 


It follows from the Lax-Wendroff theorem [9] that u N {x, t ) cannot converge boundedly a.e. 
to any function u(x, £), for otherwise u = u(x,t) would be an entropy conservative weak 
solution of (2.4), satisfying the entropy equality 


( 2 . 11 ) 


l=0; 
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this is incompatible with the entropy condition (1.1b), once the initial data u(x, 0) = 
limjv-»oo «/v(x, 0) permit discontinuous solution u{x, t). Moreover, the bounded solutions 
of (2.8) cannot even converge weakly to the entropy solution, u(x, t), of (2.4). Otherwise, 
j^FD^ x e UN ^ - and hence by (l.22c) also iJ) N e UN ( x,t ) — should tend weakly to 
which in view of (2.10) leads to the same contradiction we had in the ^-dospectral Coun- 
terexample 2.2. 


Counterexample 2.4: [18]. The scalar finite-element approximation (1.24b) of (1.1a), 
induced by the quadratic entropy, Z7*(u) = |u 2 , reads 


( 2 ' 12a ) ^u N {x Vi t) + = 0, K+i = l^ =1 _f( u u+^)) d ^ 

where u„ + i (£) abbreviates 


(2.126) 


«,/+*•(£) 2 "h Up{ (x^i, t)] + £[ujy(x,/.).i, t) u^(x,/,t)]. 


Multiplying this by ujv(x„,t) we obtain after rearrangement 
(2.13a) 

2dft UN ^ Xv ' t ^ i + Ax Ujv ^ X|/,< ^ /l,/+ 5' ~~ S 2dt^ UN ^ X>/,t ^ 2 + Ax^+^ ~ = °> 


where the consistent entropy flux F* +l is given by [18] 


(2.136) F* + i = -\u N [x v ,t) +u N (x l/+1 ,t)}h* t/+ ^- -[J f(u)du + J f(u)du]. 

In particular, this implies a global quadratic entropy conservation in time 

•2*- 2N aw 1 

(2.14) / U(u N (x,t))dx = ^2 U(u N (x„,t))Ax = ^ U(u N (x v , 0))Ax, U(u) — -u 2 . 

!/=0 >>— n 2 


u=0 


Hence u^(x, t) cannot converge boundedly a.e. to any function u(x, t), which otherwise, by 
the Lax-Wendroff theorem, would be an entropy conservative solution of (l.la). Moreover, 
by the same argument as before, bounded solutions, un(x, t), cannot even converge weakly 
to any conservative solution of the GNL equation (l.la), for otherwise the convergence 
should be strong-a contradiction. 


The essential ingredient behind the failure of convergence demonstrated in counterex- 
amples 2. 1-2.4, is the lack of entropy dissipation . Namely, in each case we have found 
an entropy function, U(u), such that the total amount of entropy (u N (x u ,t)) Ax is 

conserved in time 

2 N 2 N 

(2.15) ^2 U{ujf{x v ,t))Ax = ^2 U(u N (x „,0))Ax. 

u=0 u=0 
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This rules out (even weak) convergence to the entropy solution of our problem. In fact, in 
the finite-difference and finite-element cases, we found a local (cellwise) entropy conserva- 
tion, consult (2.9), (2.13), which prevents the existence of any strong limit as well. 

We conclude that some sort of global entropy dissipation is necessary for an consis- 
tent discrete approximation of (1.1a) to converge, for otherwise, the lack of such entropy 
dissipation is inconsistent with the augmenting entropy condition (1.1b). 
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3. CONSISTENCY WITH THE ENTROPY CONDITION 


Let U (u) be a strictly convex entropy function (we assume that the system (l.la) is 
equipped with at least one such entropy function). Multiplying (1.4) by U'(uN(x,t)) we 
have 

d d d 

(3.1) — [£7(^(2, f))] + —[F(ii N (z,t))] =< U'(u N (x,t)), —E f N u N (x,t) > . 

A possible attempt to define consistency with the entropy condition would be to require 
that the righthand-side of (3.1) is nonpositive for all IPs associated with (l.la), in agree- 
ment with (1.1b). This means that we prohibit any local entropy production by our 
approximation. In this context one is led to the concept of E-schemes [13], [17], [15]. We 
recall that in the scalar case, E-schemes are convergent in view of their Total- Variation- 
Diminishing (TVD) property [13], [17], [15], and that in the GNL case of strictly hyperbolic 
2x2 systems, E-schemes are convergent provided their solutions remain with sufficiently 
small variation [1, p.33]. Unfortunately, such an E-consistency requirement which prohibits 
any local entropy producing waves, restricts our approximation to first-order accuracy [13], 
[15]. As noted by DiPerna, however, one can allow entropy producing waves to be mixed 
with entropy dissipating waves and still retain convergence, as long as the sum of the total 
amount of entropy production, ||Up ro d(ujv(2:» t))||i,i (t), where 

r 2 jt g 

(3.2 a) Uprod(tMz)) = ~ J < ^C7'(iUtf(x)),E^tytf (a) >” dx > 

plus the total amount of entropy dissipation, ]|t7di ss (u.v(z, t)\\ L i (q, where 

(3.26) U diss (w N (x)) = < ^-U'(w N [x)),E f N w N (x) > + dx , 

remains uniformly bounded w.r.t. N. To this end, let us pick one strictly convex entropy 
function, say t/*(u), and integrate by parts (3.1) with U[u) = U*(u), obtaining 

(3.3) U,(u N (x,t))dx+ (j^Ul(u N (x,t)),E{fU N (x,t)\ = ° . 

The second term on the left is the amount of entropy dissipation rate. The counterexamples 
given in the previous section show that in order to agree with the entropy condition (1.1b), 
this amount of entropy dissipation rate should stay larger than a sufficiently small lower 
bound. In our next definition we precisely quantify the size of such lower bound. We make 

Definition 3.1: ( Consistency with the entropy condition) . The approximation (1.3) 
is consistent with the entropy condition (1.1b) w.r.t. the “relevant” class U , if for all xvw(x) 
in <f>N which are uniformly bounded w.r.t. N, there exist 
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(i) at least one “preferred” entropy function, U t {u), 

(ii) positive sequences, j 0 and 0 < r) N < rj^, 

(iii) a nonempty “relevant” class, U, of strictly convex entropy functions, 
such that for all Z7’s in U the following estimate holds 

(3.4) — ||^u;^(x)|| 2 + -^£/p rod (uMaO) < (S-K{^N{x)),E f N w N (x) \ + ||t^jvM|| 2 . 

e N r) N \ox J 

Remark (on the entropy consistency definition): The essence of the Definition 3.1 lies 
in the so called “relevant” class U. As we shall see later on, consistency in the sense of 
Definition 3.1 +L°° stability implies convergence to a limit solution of scalar and some 2x2 
approximations (3.1). This limit solution agrees with the entropy inequalities (1.1b), pre- 
cisely for those entropy pairs (U(u),F(u)) with U(u) belonging to U. Thus, Definition 3.1 
yields consistency w.r.t. the entropy functions in the class U, which motivates its descrip- 
tion as the “relevant” class. Clearly, the larger the size of U is, the sharper convergence 
results can be deduced. 

An immediate consequence from this kind of entropy consistency definition is 

Lemma 3.2: If the approximation (1-8) is consistent with the entropy condition (l.lb), 
then it is also consistent with the conservation law (1.1a). 

Proof: Using the Cauchy-Schwartz inequality to upper bound the first term on the 
right of (3.4), we obtain 

(3.5) i-|| < £7” • „,„(*) II* + ||tu w (x)f , 

with = ||17 ,, (iyjv(a:))||£,oo( J .). Multiplying (3.5) by 2e N we find 

(3.6) WE^wn^W 2 < t/" • e 2 N \\~w N {x)\\ 2 + 2e N \\w N {x)\\ 2 . 

Finally, thanks to conservation we can restrict our attention to w^[x) with zero mean, 
A- fo* ww(x)dx = 0, for which we have the Poincare inequality 

(3.7) ||t^(x)H 2 < ||~iMz)!| 2 . 

The consistency estimate (1.7) with c N = ^U^e 2 N + 2e N — » 0 follows from (3.6) and 

(3.7) . □ 


14 



Thus, Definition 3.1 implies consistency with both parts of the initial-value problem 
(1.1a, (1.1b) - the conservation law (1.1a) and its augmenting entropy condition (1.1b). 
We shall therefore refer to approximations which fulfill Definition 3.1 as consistent with the 
initial-value problem (1.1a), (1.1b), or simply as consistent approximations . 

We observe that our consistency requirement (3.4) places a rather weak restriction 
on the approximation (1.3), in that it allows for a mixture of entropy producing and 
entropy dissipating waves. Instead, it aligns with the necessary condition for convergence 
discussed already in Section 2, which requires some sort of global entropy dissipation. 
Indeed, the total mass of entropy production plus the total mass of entropy dissipation 
remains uniformly bounded in this case, as told by the essential 

Lemma 3.3: Consider the approximation (1.3) which is consistent with the initial- 
value problem (1.1a), (l.lb), and assume that it is L°°-stable, i.e., that there exists a 
constant M, (independent of N ), such that 

( 3 - 8 ) \\uN{x,t)\\ L r oc lx,t) < Moo- 

Then there exist constants (independent of N), such that 

I. The local error, E^UN(x,t), satisfies 

(3.9) ||££u N (x,t)|| i *j Iit) < Consti • s/eff -> 0. 

II. For all entropy functions U(u) in U, we have 

f T 

(3.10) j ^U pTod (u N (x,t))dt < Const n • p N . 

III. For all entropy functions U(u) in U, we have 

(3.11) || < — U'{u N {x,t),E f N u N {x,t) > \\ L l c (z,t) < Constui . 


Proof: Applying the entropy consistency estimate (3.4) to ty^(-) = u N (',t), we find 
after temporal integration that for any U (u) in U the following estimate holds 


(3.14) 


e N ft =0 \\E f N UN{x, t)|| 2 dt + ^ f l=0 U pro d(uN{x,t))dt < 

^ /<= o [h Ul{u N (x,t)),E f N u N {x,t )) dt + tf =0 \\u N (x,t)\\ 2 dt . 
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To upper bound the righthand-side of (3.14), we integrate (3.1) by parts, obtaining that 
for aU entropy functions U(u), we have 

(3.15) f (■^-U'(u N (x,t)),E f N u N (x,t))dt=[ U(u N (x,t))dx <2U 00 . 

Jt = o \ax ) Jx = o t=T 

Using this with U(u ) = U t (u), implies that the righthand-side of (3.14) does not exceed 

(3.16) 217.00 + T • , U.oo = ||U»(uAr(x,t))|U“ (x,t) - 

Hence (3.9) and (3.10) follow with Const j = Const// = 2 U too + TM Finally, according 
to (3.15) we have for all entropy functions, U(u), 

(3.17) 

r T rT rT ( d f 

J U di33 (u N (x, t))dt-J ' U prod (u^(z, t))dt = J ( ( — U'(u N {x, t)),E J N u N {x, t ) 

and hence (3.11) follows with Const/// = 2[f7oo + Const// • 7700 ], 

|| ^ U (tl//(x, tj), Ej/Uff(x, t) > = /(=o[^diss( u i\r(x, t)) + Up rod( u w(x, t))]<ft < 

^ 2IJ/OJ, + /f = o Up ro( j(u^/’(x, t))df] ^ C onst/// . 

Remarks : 

1. The estimate (3.9) shows that the local error of consistent schemes, E^u N (x,t), 
tends to zero independently whether the underlying solutions are smooth are not. 

2. The estimate (3.10) shows that for any entropy function U(u ) in the “relevant” class 
U, the total mass of its production remains bounded (by Const//'T} 00 ). Moreover, the total 
mass of its entropy production tends to zero in case 

r T 

(3.19) / U prod (u N (x, t))dt < Const// • r] N — >• 0 . 

Jt= o 

3. According to the first remark 

(3.20a) F N (u N (x,t )) = i^u^x,*))- < U\u N {x,t)),E f N u N {x,t) > , 

is a numerical entropy flux which is consistent with the differential one 

(3.206) ||F(u*(x,i)) -F>(u;v(x,t))|| < ||U'(u 7v (x,t))|| £ « (M) • \\E f N u N {x,t)\\ — ► 0 . 

The estimate (3.11) shows that for any entropy function U(u ) in the “relevant” class U, 
the total mass of its production + dissipation, || ~^[U(un(x, t)] + f^\F//(u//(x, £))]|U, l oe (* I t) 


) dt < 2U oa , 
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remains bounded, as asserted. This agrees with DiPerna observation [l, p.39] that the 
control on the total entropy production is independent of the choice of the preferred entropy 
function 77* (u). 

4. Let us restrict our attention to consistent approximations of the scalar initial-value 
problem (1.1a), (l.lb), and apply Gronwall’s inequality to (3.3) with 7/*(u) = |u 2 , 

1 d ( d t \ r2x 

(3.21) u 2 N (x,t)dx = - I —u N (x,t),E f N u N (x,t) ] < u%(x,t)dx . 

L dt J z=0 \^C/X J * z=0 

We conclude that independently of the Instability assumption (3.8), we have 

(3.22) ||uAr(M)|U*(x) < Af 2 • 

In case 77* (u) is a strictly convex entropy function, the same conclusion also holds for 
systems of conservation laws, by appealing to (3.10) with the quadratic part of 77* (u). 

The intricate point in our entropy consistency Definition 3.1, is the size of the “relevant” 
class U. The optimal situation of course, occurs when U = l/ a ii which includes aU the 
admissible strictly convex functions associated with (l.la). In this respect, it would be 
useful to have a sufficient consistency criterion which guarantees the optimality of U by 
checking only one preferred entropy function, say U*(u). Our next result provides us with 
such a convenient criterion. 

Lemma 3.4: ( Sufficient criterion for consistency) . Assume that for all wn[x) in <f>n which 
are uniformly bounded w.r.t. N, there exist (at least) one strictly convex entropy function, 
17* (u), and vanishing sequences, c^ I 0 and dx j. 0, such that the following two estimates 
hold: 

A. ( Consistency with the conservation law) 

Q 

(3.23a) li^ u ’iv( x )l| < CNll^tM*)!! > 

B. ( “Enough” entropy dissipation) 

(3.236) d N \\^-w N {x)\\ 2 < ^-Ul(w N (x)),E(fW N {x)^J + ||u; J v(x)|| 2 . 

Then the approximation (1.3) is consistent w.r.t. all admissible entropy functions, U = U a \\, 
provided 

(3.24) < Const. 
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Proof: By Cauchy-Schwartz inequality and (3.23a) we have for aH entropy functions, 
U(u), 


(3.25) U prod (w N (x)) < K • c N \\-^w N {x)\\ 2 , U Z = 1 1 U"{wu ( 2 ) ) II £«= (x) . 

Take e N = 2 Const. 2 d N j 0 and rjs = < f?oo = 2Const.U^ ) ; then using (3.23a), 

(3.25) and (3.23b) we obtain for all strictly convex entropies, U(u), 


^l|£w t M*)|| 2 + ^tfprodKrfc)) < 2Con]tld N C N 


h w **( x ) 


+ 2 dux, U 'L C N 


h w N{x) 


< 


< $[d N + d,N 



< 


< (^Ul{w N {x)),E^w N [x)) + ||wjv(x )|| 2 . □ 


Remarks: 

1. The sufficient consistency condition provided in Lemma 3.4 is sharp, in the sense 
that if (3.23) holds, and if instead of (3.24) we have 


(3.26) 


C N 

d,N 


0, 


then approximation (1.3) does not converge. Indeed, if a limit solution of (1.3) exists in this 
case, then it is necessarily t/*-entropy conservative, for (3.23), (3.26) imply that both the 
U t -entropy production and entropy dissipation tend to zero. Put differently, approximation 
(1.3) is not L°°-stable in such case, for otherwise it would contradict convergence along 
the lines of the scalar counterexamples in Section 2. 

2. The entropy dissipation estimate, (3.23b), with no “lower-order” term ||twj V (x)|| 2 , 

(3.27) d N \\-^w N (x)\\ 2 < ^-Ul(w N [x)),E^w N (x)^ , 


can be viewed as an accretivity condition on the discrete local error E * N . However, it is 
worthwhile noting that such accretive approximations are restricted to first-order accuracy 
due to the following argument: 

3. By the weighted Cauchy-Schwartz inequality, the righthand-side of (3.27) is smaller 
than 

eC|||w*)II ■ K®*(*)ll . 


a N | 

U" 1 
^*00 


dx 


«>jv(a;)|| < ||^i»jy(*)|| • 


and hence (3.27) implies 
(3.28) 
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The inequality (3.24) implies that > Const. • jj. Consequently, the inequality (3.28) 
shows that approximations which satisfy the entropy dissipation estimate, (3.27), are re- 
stricted to first-order accuracy. Instead, we adopt here the weaker entropy consistency 
estimate, (3.4), which, roughly speaking, requires the total entropy dissipation minus pro- 
duction to be proportional to the local discrete error E^wn(x). This will enable us to 
deal with higher (than first) order approximations of (1.1a), including the modern high- 
resolution finite-difference schemes [2], [3], [14], [15]. 
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4. CONSISTENCY OF THE SCALAR VANISHING VISCOSITY METHOD 


The vanishing viscosity method is one of the most effective ways to provide affirmative 
answers to the questions of consistency and //“-stability, in the context of nonlinear sys- 
tems of conservation laws. To this end, one starts from a given basic approximation , P^, 
which is consistent with the conservation law (1.1a) 

(4.1a) ||£jtM*)|| < c N \\~w N {x)\\ 

Such basic approximation is then appended with a so called vanishing viscosity approximation , 
- enQn£z , which results in the modified viscous approximation 

(4.16) H = Pn~ 

Here e N j 0 is the vanishing viscosity amplitude, and Q N : L°°[0,27r] — > $ N , is a possibly 
nonlinear, / dependent, spatially bounded operator 

(4.2) WQnj^nW < Qoo||^^(a:)|| • 

The boundedness of Qn implies that —SnQn§^ is consistent with the zero flux, i.e., 

d d 

II - s n Q n —w n (x ) || < Q TO £iv||^Wiv(a:)|| > Qoo^n -> 0 ; 
consequently, the discrete local error of the modified viscous approximation (4.1b), 

(4.3 a) E f N w N (x) = f{w N (x)) - P f N w N (x ) = E f N w N (x) + e n Q n ~w n (x), 

satisfies, by (4.1a) and (4.2), 

Q 

(4.36) ||^Ju>iv(*)|| < Civ||^^iv(2)|| , ~c N = c N + Qco£n -> 0 . 

In this manner, the vanishing viscosity method retains the consistency with the conserva- 
tion law, (1.1a), of its underlying basic approximation. 

There are various ways to tune the viscosity parameters so as to guarantee the full 
consistency of the viscous discretization (4.1b). Namely, we seek for e n and Qn so that 
Pjy = p/f — enQnj^, which remains a consistent discretization of the conservation law 
(1.1a), gains in addition, the consistency with the entropy condition (1.1b). Taking ad- 
vantage of the scalar case, we make the canonical choice of U t (u) = \u 2 as our preferred 
entropy function for the sufficient entropy consistency condition stated in Lemma 3.4, 
which is then simplified to 
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Theorem 4.1: Let P^ be a consistent discretization of the scalar conservation law 

( 1 . 10 ) 

Q 

(4.4) \\E f N w N {x)\\ < c N \\—w N {x)\\ , c N 0 . 

Assume that there exists a vanishing sequence, d n j 0, such that for uniformly bounded 
w N (x) in 4> n , 

(4.5) d N \\^-w N (x)\\ 2 < ^-w N (x),E f N w N (x)^j + ||tu;v(x)|| 2 • 

Then we have 

(i) The vanishing viscosity method, P f N = — zQnjz, ts consistent w.r.t. all strictly 

convex entropy functions, U = U a \\, provided 


(4.6) 


- — ^ Const . , Cft = cn + Qoo^n • 

d N 


(ii) The total amount of entropy produced by the basic approximation, Pfj, 

f T f T d f 

(4.7 a) / U pmd [u N (x,t))dt~ / -<—U'(u N {x,t)),E I N u N [x,t)>~dxdt, 
Jt =0 Jt=0Jx—0 ox 

tends to zero, provided 


(4.76) 



Remark : The second part of Theorem 4.1 shows that due to the presence of vanishing 
viscosity parametrized according to (4.7b), the entropy produced by the basic approxima- 
tion tends to zero. Indeed, by (4.4), (4.5), (3.15) and (3.22) we have 


(4.8) 


- St = 0 lx =0 < $Z U '( u n{x, t)),E f N u N (x, t ) > dx dt 


U'L 


cn 


JZ u N{x,t) 


2 

< 


<P” -£(r+i)MJ — .o. 


The following consequences of Theorem 4.1 are at the heart of the matter. 

Corollary 4.2: Let Pj) be a consistent approximation of the scalar conservation law 
(1.1a), 

(4.9) ||^jvWjv(*)|| < c N \\ — tu w (a:)|| . 
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Assume that for all w^(x) in there exists a positive constant (independent of N), 
Const. > 0, such that, 

( d d \ d 

—w N (x), Qn-q^wn{x)] + lktf(a:)|| 2 > Const.e N \\— to^(s)|| 2 . 

Then, the resulting vanishing viscosity method 

# 

d d d d 

(4.11) —[u N (x, t)] + —\P f N u N {x, t)] = s n —[Q n —u n (x, t)] 

satisfies 

(i) It is a consistent approximation of the scalar initial-value problem (1.1a), (1.1b), w.r.t. 
all strictly convex entropies, U = Z/ a u, provided 

(4.12a) — < Consto < Const. . 

£N 

(ii) The total amount of entropy produced by its basic approximation in (4.7a) tends to 
zero, provided 


(4.12 b) — — ► 0 . 

Proof: By (4.9), and (4.12a) or (4.12b) with sufficiently large N, we have 

(fe w w(*), £&«;*(*)) > -||^ta^(i)|j • ||^it» JV (x)|| > -c^||^u;^(z)|| 2 

> -Const. eE N \\^w N {x)\\\ O<0=§^<1. 

Hence, (4.5) is fulfilled with d.^ = Const. (1 — 0)en, for 

dAr||^Ar(®)|| 2 = Const.e N \\j^w N {x)\\ 2 - Const.0er, v ||^i£; iv (x)|| 2 < 

(4.13) < [sn (-tw N {x),Q N £w N {x)) + ||^iv(x)|| 2 ] + [f^w N (x),E f N w N {x)) = 

= {£ w N{x),E f N w N (x)) + ||u>at(z)I | 2 . 

Now, with this choice of d/v we have 

Cn _ 1 CjV + QqqSn < 0 , Q 00 . 

~d N Const. (1 - 6)e n ~ 1 — 6 Const.( 1 - 6) 
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moreover, (4.12a) implies 


1 


cn 


cn 


d N Const. (1 - 6)£n 
and the result follows from Theorem 4.1. □ 


10 , 


The first part of Corollary 4.2 guarantees consistency, provided the amount of entropy 
produced by the basic approximation, Pjy, is dominated by the amount of entropy dissi- 
pated by = Ph — £nQn§^- In case the basic approximation was entropy conservative 
to begin with, then we can do with even less vanishing viscosity as told by 

Corollary 4.3: Let P^ be a consistent approximation of the scalar conservation law (1.1a), 
which is (quadratic) entropy conservative, i.e., 

(4.14) (^-w N (x),E f N w N (x)^j = 0 . 

Assume that (4-10) holds. Then the resulting vanishing viscosity method ( 4-11 ) is a con- 
sistent approximation of the initial-value problem (1.1a), (1.1b) w.r.t. all strictly convex 
entropies, U = U a \\, provided there exists an arbitrary positive constant (independent of N), 
Consto > 0, such that 

(4.15) — < Consto » Consto > 0 . 

£n 


Proof: Choosing d n = Const.ex j 0, then (4.14) and (4.10) yield 


(jzWN(x),E f N w N (x)) + ||u>w(z)|| 2 = e N (£ w n {x),Q n £w n (x )) + ||u> N (z)|| 2 > 


>d N \\fw N (x)f 

and the result follows from Theorem 4.1 since 


cn 

d N 


on 4 " Qoo^n 
Const.Es 


< 


Const. 


(Consto + Qoo) • 


□ 


Corollaries 4.2 and 4.3 show that the vanishing viscosity method enables to gain en- 
tropy consistency while retaining consistency with the conservation law. It is possible 
however, that this method gains the entropy consistency at the expense of lowering the 
order of accuracy of the underlying discretization P^. The vanishing viscosity method 
should be carefully parameterized in order to retain both the entropy consistency as well 
as the original order of accuracy. An extreme situation in this respect is provided by the 
“infinitely-order” accurate spectral and pseudospectral methods, which bring us to 
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The Spectral Viscosity (SV) Method [20] [10]. We consider spectral vanishing viscosity 
modifications, defined in terms of a convolution with a 7r^-kernel, Q^(2:), 


d d 

(4.16a) -£ N Q N —w N (x) = -£ N Q N {x) * — Wn ( x ), 

where Qn{x) = Qn( w n(x), x) is a possibly nonlinear viscosity kernel of the form 

(4.166) Q n {x) = £ Qkt ik \ 0 <Q k <Qoo. 

k=-N 


Appending this either to the spectral ( P n = S^) or ^>dospectral (P# = iJjn) approximation 
of (1.1a), results in the SV method 

d d d d 

(4.17) — [u N [x, if)] + — [Pn/{un[x, t))] = £n-^[Qn(x, t) * — uw(x, *)]; 

it can be efficiently implemented in the Fourier space as, see (l.ll) or (1.13), 

(4.18) ^ujt(f) + = ~ £ N^ 2 Qk(t)u k {t). 

According to Theorems 7.1 and 7.3, the spectral and pseudospectral methods are first- 
order accurate with (1.1a), i.e., they satisfy (4.1a) with c ^ = Const 0 • i, where Const 0 = 
(* + js) II f'(u N {x, f))l|i,- (*,«)• Appealing to Corollaries 4.2 and 4.3 we arrive at 

Theorem 4.4: ( Consistency of the SV Methods) The following SV methods are con- 
sistent approximations of the initial-value problem (1.1a), (l-lb) w.r.t. all strictly convex 
entropies U(u): 


I. The spectral viscosity approximation 

d d d d 

(4.19a) — [u N (x,f)] + ^[SV/(u;v(x,f))] = s N —[Q N (x,t) * —u N (x,t)}, 

with spectral viscosity parameters satisfying 


(4.19 b) 



Const. — 


1 



— <e N 10, 


Const. > 0 . 


II. The ipdospectral viscosity approximation 

d d d d 

(4.20 a) —[u N {x, t)] + -7^[ip N f[uN(x, 0)] = £n^;[Qn{x, t ) * —u N (x, f)], 

with spectral viscosity parameters satisfying (here Const 0 = (l + \\f'\\L~ c {x,t)) 

(4.206) Q k (t) > (Const. - < s n I 0 , 


Const. > Consto . 
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III. The total amount of entropy produced by the basic spectral approximation (where 
Pn = Sn), or the ipdospectral one (where Pn = )> tends to zero, 

(4.21 a) - f f < ^-U'(u N (x,t)),[I - P N ]f(u N (x,t)) >~ dxdt — > 0 , 

Jt=0Jx=0 OX 

with spectral viscosity parameters satisfying 

(4.2U) <?* W >(con 5t .-^) + , ^-*0. 

Proof: According to Corollaries 4.2 and 4.3, consistency follows if (4.10) holds, and 
by Fourier transform this is equivalent with 

(4.22) £Nk 2 Qk{t) + 1 > Constant: 2 . 

In the spectral case there is no entropy production 

(j^w N (x),[I - S N ]f(w N (x))^ =0 

and the first part of the theorem follows from Corollary 4.3. In the ^dospectral case we 

V 

apply Corollary 4.2(i) and the second part of the theorem follows. Finally, the third part 
of the theorem follows from Corollary 4.2(ii). □ 

Example 4.5: The following choice of spectral viscosity for the SV method (4.17), 

d d i N wP 

(4.23a) e N —[Q N * —u N {x,t)\ = £ k 2 Q k {t)u k (t)e' kx , m 2 N = , 

|fc|=m w ''' 

where 

- N& 

(4.236) Qjfc(t) > Const. — -p-, \k\ > m n , 0 < (3 < 1 , Const. > Const 0 , 

demonstrates a consistent approximation of (1.1a), (1.1b) w.r.t. all strictly convex en- 
tropies, which retains the “infinite-order” accuracy of the underlying spectral or pseu- 
dospectral method. Moreover, if 0 < (3 < 1, then the entropy produced by the basic spec- 
tral or t/>dospectral approximation tends to zero, when augmented with (4.23a), (4.23b). 
For example, the SV approximation 

d d 1 N (k — \/N) 2 

—[u N {x,t)} + —[P N f{u N {x,t))} = -Const.— k 2 -——y=-u k {t), 

|Jfe|=v/W iV ‘ v/iV 
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is ‘infinitely-order’ accurate approximation of (1.1a), and at the same time it is consistent 
with (1.1a), (l.lb). 

Remark : Let us rewrite the SV method (4.17) in the form 
(4.25) 

d d d d d 

— [u N (x,t)\ + — [f(u N {x,t))]-—[(I-P N )f(u N (x,t))} + £ N ~ Q N Q^u N (x,t) . 

This highlights the fact that the S V approximation has two sources of errors: the first term 
on the right is the error committed due to spectral discretization; the second term is due to 
the presence of spectral viscosity. The viscosity parameterization in (4.19b), (4.20b) guar- 
antees that the local discretization error is dominated by the dissipative spectral viscosity 
error. 

Next, we turn to study the vanishing viscosity method in the context of 


Finite-Difference and Finite-Element Methods . We consider conservative approxima- 
tions of (1.1a), whose numerical flux ^ +1) is given by (abbreviating 
Au v+ l(i) EE u N (x v+1 ,t) - u^(x^,t)) 

(4.26) h f v+l _ = f{u N {x u ,t )) + f[u N {x v+1 ,t))\ - ^Q v+ i_Au u+ i(t) . 

The conservative approximation (1.14) then recast into the form 


£ t u N (x,,,t) + 2 ^[/(ujv(a: J /+i,*)) - /(«jv(*i/-i,<))] = 

= — Q„_iAu„_i($)] , 


which reveals the role Q v+ 1 play as the numerical viscosity coefficients in (4.27). Simple 
linear examples are provided by centered finite-difference method (1.19), where 


(4.28) Q v+ i = 0 , 

and the piecewise-linear finite element method (2.12), whose numerical flux 

(4.29a) ^„ + i = ^ 0)] — ’ 

corresponds to the viscosity coefficient, Q u+ i = Q* i, consult [18], [21] 

> 2 2 

(4.294) <?* +J = 2 - J^_ l/2 (j - ?) • 
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Once more, we would like to emphasize that our present treatment of the vanishing viscosity 
method allows for nonlinear viscosity coefficients, Q u+ 1 , which may depend on the com- 
puted solution Q v+ % = Q„ + i(x,f) = Qv+i{uN(x,t),x,t] f). I n particular, this includes the 
essentially nonlinear features which characterize the modern nonoscillatory finite-difference 
schemes [3]. 

In order to interpret (4.27) within the vanishing viscosity framework (4.1a), (4.1b), 
the following viscosity operator Q N : L°°[0,27r] —* n N is introduced: We construct the 
fl-jv-polynomial 

K N (x) = K N (x-,Q) = £ K k e ikx , 

I brs-JV 

which interpolates the Q-weighted cell averages, tF(x) = — / i _ £ w(£)d£ of tu(x), 


2 N 


(4.30a) K n {x = x„ + i) = Q v+ lw[x v+ i) , K k = ^ £ Q„ + »tF(x„ + »)e tk *"+h , 

Z7r v=0 

and then define Q^w(x) as the sliding average of this interpolant 


(4.306) 


_ _ 1 /•!+ — 

Q n w{x) = K n {x\ Q) , K n {x) = — y _ J K n ($ ; Q)d£ . 


By applying such averaging plus summation by parts to JC^(x,t) we find, in analogy with 
(1.17), 


(4.31a) 

where 


N 


Q N u N (x,t) — ^2 qk{t)e tkx , 

k=-N 


■i 2 JV 

(4.31.6) ft(0 = 2 ^ - Q v _iu JV (x^_.,t)]e-’ 7:s> ' . 


Now, since (x = x„ + i) = u„ + i(t), we have 


(4.32) 


d 1 2N 1 


27 rffc Ax 


and we conclude that |A xf^QN-§^UN(x,t) is the ^-polynomial which interpolates at the 
equidistant point x„ = i/Ax, the viscous part of (4.27), 


1 q d 1 

(4.33) 2^ X ^x^ N ^x U n ^ X = = 2Ax^ |,+ ^ U|/+ 5'^ _ ’ 
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Thus, the conservative approximation (4.27) can be viewed as the vanishing viscosity 
method (4.1), which is based on the consistent centered finite-difference approximation 
(1.19) and augmented by the vanishing viscosity approximation (4.33), 

(4.34) Pfj = FD Ax f(-) — e N Q N — , e n = . 


The representation of a given approximation as a vanishing viscosity method is in general 
not unique: one can assign any part of it to serve as a basic approximation (as long as it is 
consistent with the conservation law), and then consider the rest as the vanishing viscosity 
contribution. For example, we can represent (4.27) as a vanishing viscosity method which 
is based on the finite-element approximation, P ^ = FE Ax f(-) in (1.26). To this end we 
use (4.29a) to rewrite (4.27) in the equivalent form 


(4.35a) 

where 


dt u N(xv,t) + A(r [h„ + i h u _i\ - 9AT [- c, i/+|Att 1/+ |(t) D v _iAu u _^(t)] 


Ax 


2Ax 


(■*•356) D„ h = Q r+i - <?- + . , 

are the numerical dissipation coefficients of (4.27). Now, by introducing the corresponding 
viscosity operator D n ' 2/°°[0, 27t] — > ir N , 

(4.36) D n w(x) = K n (x\ D) , K n (x = x v +^\D) = D u+ vw{x vJr P) , 

we conclude that (4.35) is associated with the (modified) discretization operator 

(4.37) ~P f N = FE Ax f(-) - e N D N ~ , e N = ~. 

Which representation should we prefer? Corollary 4.3 suggests that in order to obtain 
sharp consistency estimates for a given approximation, the key step lies in isolating an 
entropy conservative part of it as a basic approximation, and considering the rest as the 
vanishing viscosity contribution. We therefore prefer (4.37) over (4.34) to make 


Theorem 4.5: The finite-difference/element approximation 


(4.38a) 


ft U n[x u , t ) +23^[/(uiv(x 1/+1 , t)) - f(u N (x v -!, *))] = 

= jkzlQ V+ 5^ U v-\- ~ < ?i/-|-^ U v-i(0] > 
with viscosity coefficients, Q u+ i, satisfying 


(4.386) 0 < Q* v+ i + Const. < Q u+ i_ < Qoo , Const. > 0 , 
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is a consistent approximation of the scalar initial-value problem (1.1a), (1.1b) w.r.t. all 
strictly convex entropies. 

Proof: We first verify the (quadratic) entropy conservation of the basic finite-element ap- 
proximation, Ph = FE Ax f(-). By the exactness of the trapezoidal rule for ^^-polynomials 
and the identity (2.13a), we have 


(»«Wi §;FE^ z f(w N (x))) = Ejjfo “»(*••) - K.-il ' Ax = 


= ESoK.1 - K- i] = 0 , 


and since < lOjv(^c) j Jif{ w N{ x )) >= j^F(wn(x)) is a perfect derivative, 


^w N (x), -^-f(w N (x))^j = 0 . 

Integrating by parts the difference between the last two equalities, we conclude 

(4.39) ^~w N (x),E f N w N {x)j = (J^ w n{ x ), [I - FE Ax \f(w N (x))^ = 0 . 

This reaffirms the quadratic entropy conservation, (4.14), of the basic finite-element ap- 
proximation (2.12), which was already indicated (with different terminology) in Coun- 
terexample 2.4. According to Theorem 7.5, the finite-difference approximation (4.38) is 
first-order accurate with (1.1a), i.e., (4.1a) is satisfied with c n = Const 0 • -h, and hence 
(4.15) holds. 

Appealing to Corollary 4.3, consistency is therefore guaranteed if there exists a positive 
constant such that the following estimate holds 

— ( g^w N (x),D N —w N (x) \ + ||wn(s)|| 2 > -Const.Ax\\—w N (x)\\ 2 , Const. > 0 ; 
by the trapezoidal rule this is the same as 



and the result follows from our assumption (4.38b), for 


D v+ 1 = Q u+ i - Q* + ± > Const. > 0 . 
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Remarks : 

1. We observe that the total amount of quadratic entropy dissipated by (4.38a) equals, 
in view of (4.39), 


= X ( jk W »( X )- D H'k W H ( X )) = 

(4.41) . 

which justifies calling D v+ ±. the numerical dissipation coefficients of (4.38a); consult [18] 
[21] for the general case. 

2. The finite-difference/element approximations (4.38a), (4.38b), are restricted to first- 
order accuracy, in agreement with Corollary 7.6. Indeed, the proof of Theorem 5.4 makes 
no use of the “lower-order” term ||iu;v(:e)|| 2 in (4.40), hence (3.27) holds. 

Theorem 4.5 tells us that the conservative approximations (4.27) are consistent with 
the initial-value problem (1.1a), (l.lb), provided they contain uniformly more viscosity 
than the entropy conservative ones, (1.24), (1.25). To obtain a sharper localized version of 
such result, we shall prefer to represent these approximations in terms of local ‘hat’ basis 
functions, {A*.(x)}^ = _ // , instead of the global trigonometric ones, {e , * z }£L_ JV , we have used 
so far. To this end we may proceed as follows. 

We first recall that in (4.35), the finite-difference/element methods (4.27) were based 
on the finite-element approximation with mass lumping (1.24b). Before we turn to consider 
the viscosity contribution in (4.35), however, we shall prefer to concentrate on a slightly 
different, yet closely related basic approximation - the finite-element approximation with 
no mass lumping in (1.24a). Example 1.7 provides us with the piecewise-linear formulation 
of this basic approximation: it is associated with the piecewise parabolic discretization 
operator in (1.35), P^w N (x) = FE Ax f[w N [x)), such that 

(4.42) ^-P^w N {x) = £ ^[ C '“ lA 7l* A *( a: )- 

The important ingredient of this approximation, is the orthogonality of its discrete trun- 
cation error, ■^Ej V wn(x), to the ^-space, 

(£*&»»(*). M*)) = (£/(«w(*)).M*)} - &[c-‘A7J* (Afc(z),A„(x)) = 

= 17,+ 1 - 7,-jl - £[mc-‘a 7), = o . 

This brings us to the canonical formulation of the basic finite-element approximation (1.36) 
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as a Galerkin method, namely, u^(x,t) is the piecewise-linear approximant such that of 
all <Pn(x) in <f>N we have 

(4.43) (vn{x), u N (x,t ) + ^[/(«jv(x,t))]^ = 0, <Pn(x) <E 4> N . 

Counterexample 2.4 shows that this method is inconsistent with the entropy condition 
(1.1b), precisely because it is quadratic entropy conservative: indeed (4.43) implies that 
(4.14) holds, for 

Wjv(x), ■^E f N w N (x)^j = 0 , w N (x) e <f>N . 

We therefore appeal to the vanishing viscosity method. In this context, we seek a piecewise- 
linear approximant, Uff(x,t) — UN[xk,t)A.k{x), such that for all <pn[ x ) in 4>n we 

have 

( d d \ A t. ( d d \ 

< fiN{x),—u N {x,t) + ^[/(Mx,*))]J + — [g^'PN{x),g^u N {x,t)J = 0. 

Here D = {-D v+ i > 0} are given numerical dissipation coefficients, so that the vanishing 
viscosity contribution enters through the weighted inner product on the left 

(v(x),u;(x)) u = £ *«V+i / < v(x),tu(x) > dx , w = {w„ + i}jL_w • 

v =-JV ■'*•» 2 

Choosing <Pn[x) — A„(x), then (4.45) reads at the collocation points x„ = uAx, 


( d ' 

—w n {x),E i n w n {x) 


(4.46a) 


dt [e^N^iz-ls t) gU^(x i/,t) + 6 u JV (x v+ x,t)] + 

"i"2 0) — f{ u N{ x t/-U 0)1 = — > 


where in view of (4.29a) , (4.29b) we have, in agreement with (4.35) , 

( 4 - 46& ) Qv+\ = Qt+i + Ah-> • 

This approximation corresponds to (the unlumped mass version of) the finite-difference/ 
element approximation (4.27), or equivalently (4.35). It is based on the piecewise-parabolic 
finite-element approximation in (1.35), 

(4.47a) P' N w(x) = £) / l+ .fWz) , A/ = (T‘A7, 

k= — N 
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which is appended by the piecewise parabolic viscosity approximation, e^ = 

where 


(4.476) 


N 


D n w(x) = Yl PHi+lfWO*) > = D k+ iw k+ i 


k=—N 


The difficulty of working with the explicit representation of the modified discretization 
operator, P f N = P* N — ^f-D N gg, lies in the nonlocal inverse of the mass matrix, C -1 , which 
enters (4.47a), (4.47b). It would be more convenient, therefore, to deal directly with the 
original weak formulation in (4.45): it asserts that for all C °°- test functions tp(x) and 
arbitrary ^-approximants, ^>jv(x), we have 


(^(a;)>ft“w(;M) + &[/(«*(*,*)]) = {<P[x) ~<Pn{x), §- t u N (x,t)) + 

(4.48) + ( [<p{x ) - <p N {x), £[f[u N {x,t )) ]) - 

-T-{& ( Pff{x),£;U N (x,t)) D . 

The novelty of the weak formulation (4.48) is that it allows us to assign to <p(x) different 

<j> Ar-approximants, (pn{x). Different choices of <Pn(z) amount to formulation of (4.46) 
within different setups. 

For example, if we take <Pn(x) = <Pn{x) to be the -^-projection of <p(x), 

N l r 

(4.49) <Pn(x) = ]C £>iAjfc(x) , (p k = C 1 • — / <p{x)h. k {x)dx , 

k=-N A. X J 

then we recover the piecewise-parabolic discretization operators in (4.47). Indeed, since 
tp(x) — <Pn(x) is orthogonal to the <£w-space, we are left with the second and third terms 
on the right of (4.48). The second term is a weak formulation for the local truncation error 
of the basic approximation 

(4.5°a) ^-E f N u N (x,t)[<p\ = (^p(x) - <p(x),-^-[f(u N {x,t))]j , 

so that Eff = [I — FE^ x \f{’), in agreement with (4.47a). The third term represents (in a 
weak form) the truncation error due to the presence of vanishing viscosity 

Ax d d Ax ( d d \ 

t 4 - 506 ) -Tai DN a; nN[x ' m = T ) D • 

in agreement with (4.47b) . The difficulty with this representation lies again in the nonlocal 
inverse, C' -1 , which enters in the definition of £>jv(x) in (4.49). 
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Instead , let us assign ^(z), as was done in [4], to be the piecewise-linear interpolant 
of <p{x), 

N 

(4.51) <Pn{x) = ^2 <p[x k )A. k (x) . 

k=-N 

We now consider the three terms on the right of (4.48). The third term is a weak repre- 
sentation for the truncation error due to the presence of vanishing viscosity 

• . \ __ d , .. , Az ( d . . d , Az 

(4.52a) e N —D N —w N {x,t)[<p\ = — — \^—<p N {x),—w N [x,t) j , e N = — . 

This corresponds to the piecewise-constant viscosity operator (here x k +i{ x ) denotes the 
characteristic function of the [z*, zjt + i]-cell) 

N 

(4.52 6) D N w N {x,t) = Y, D k+ i^ k+ i{t)x k+ L{x) , %>(<) = W(*fc+J,0 • 

k=—N 

The second term is related to the discretization of the spatial flux, f x , 

(4.5 3 a) ■^-E f N w N {x,t)\<p}= ^p(x) -ip N (x),^-{f(w N {x,t))]^j , 

where, as usual, 

(4.536) Ehw[x,t) = /(io(z, f)) — P^w[x,t) . 


This uniquely determine a piecewise-constant spatial discretization operator 
(4.53c) Pl,w{x,t) = £ 7fc+i(<)x*+i(z) , 7 k+ dt) = 7{w{x k+ i.t)) . 


We observe that -^P^w(x,t) and eNj^D^f^w(x,t) do not lie in the <}> N -space] they belong 
to the space of measures W -1,0 °. To balance these W -1, °° terms we have to “discretize” 
(in space) the temporal flux, u u as well. The first term on the right of (4.48) gives us the 
weak representation of the truncation error due to this “temporal” discretization 

(4.54a) ^E^w N (x,t)[<p]= ^p{x) -ip N (x),^w N {x,t)j , 

or equivalently, 


(4.546) E%w{x,t)[p] = {<p{x) - <p N (x),w(x,t)) . 

The discrete local error of our approximation consists now of two different sources of errors, 
(4.55a) E N Wff(x,t ) = E$Wn(x, t) + E$wu{x,t) , 
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where E$WN(x,t) is associated with discretization of the temporal flux, 

(4.556) E$w(x, t) = Epfw(x,t) , 

and E$WN(x,t) is associated with the spatial discretization (basic + viscosity approxi- 
mations) 


(4.55c) 


E$w(x,t) = E f N w(x,t) + ^—DN~-w{x,t) . 


Ax „ d 


2 dx 

The (quadratic) entropy consistency estimate corresponds to (3.4) now reads 


t N II 0 + a x ^W W N{.X,t)^\ H -l^ x ^ + ri N ^PToa[ w n{ x , 0) — 


j9 f( x ) 


(4.56 a) 


< (§i™N{x, t) , E$ w n (x, t)) + (£w n (x, t),E$w N (x, tfj + j|u;^(x,t )|| 2 . 

Here, the quadratic entropy production, U p ro d(u>w(x,t)), takes into account the additional 
discretization of the temporal flux by modifying (3.2a), 

c^prod (u>jv(x,t)) = A- 1 1 < W/^(x,t), gjEfJww(x,t) > HjJj^( Zl i) + 

(4.566) 

+11 < £wN(x,t),E$w N (x,t) >~ || ifoc(Zit) , e N | 0 . 

We now arrive at the localized version of Theorem 4.5. 


Theorem 4.6: (Upwind Differencing) . The finite-difference/element approximation 

ft [k u N(x v -i ,t) + \u N {x u ,t) + |tt^(2y+l>0] + 

+ 2 k[f( u N(xt, + i,t)) - f(u N {x v . u t)\ = 2 Zz\Q v +%&u v+ i{t) - Q U _1 Au,.i(f)l , 

is consistent with scalar initial-value problem (1.1a), (1.1b) w.r.t. the quadratic entropy 
U = {|u 2 }, if its viscosity coefficients, Q u+ l, satisfy 

[ X rxu+i l 1 / 2 

Ax /i 1/ { u n{x, t))| dx 


(4.57a) 


Proof: The proof consists of three steps. 

(i) Entropy dissipation. The total amount of (quadratic) entropy dissipation minus entropy 
production of (4.57a) equals - in view of (4.52a), (4.53a), (4.54a) and the piecewise-linearity 
of w N (x,t ) 
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(4.58) 


^w N (x, t),E$w N (x, i)) + (£w A r(x,t),E$wx(x,t)') = 

= -§- t E^w N (x,t)[<p = iUAr(x,t)] - §^E f N w N (x,t)[tp = tow(x,t)]+ 

+ = U>N(x,t)} = 

= T-{-k w N{x,t),^w N [x,t)) D = lEfL-N^+il^+iWI 2 . 
in agreement with (4.41). 

(ii) Entropy production . The viscosity approximation does not produce entropy, for 

Ax d d 1 ^ 

< -j-^vj N [x,t),D N —w N {x,t) >= - Y, D v+ i|Au;^ + i(t)| 2 x v+ i(x) > 0 . 

Consequently, the spatial contribution to the entropy production in (4.56b) is upper 
bounded by 

(4.59a) 

II < £w*(x,t),.E^W(x,t) >" || L i (s) < || < £w N (x,t),E f N w N (x,t) > || Ll(l) < 

< II &«>*(*» Olli + ^-||^u;jv(*, t)||i,-i • 

Using the “super-approximation” estimate [4, Lemma 2.1], we can upper bound the tem- 
poral contribution to the entropy production in (4.56b) 

Aril < WN{x,t),§- t E^w N [x,t) > < 

(4.596) < Ar SU P I - (wN(x,t)<p(x)) N , £w N (x,t)) | 2 < 

IMI =1 

< r?oo^f^||£u;jv(x,t)|| 2 , 7)00 = Const. 

(iii) Discrete error . The spatial contribution does not exceed 

^\\£zE$ W N{ x >t)\\ 2 H-i( x ) < 

(4.60a) < ||^i>»E w »( x ' , )lli>-i + 7 ^II £ A' ui k( i . 1 )IIi3-i < 

An upper bound for the temporal contribution is given by 
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( 4 . 606 ) 


sup \(<p(x) -<p N {x),-§- t w N {x,t))\ 2 < 

II 11=1 

< . 

Now, let us choose = Const. Ax J. 0 and t)n = Const, with sufficiently large Const. > 

0. In view of (4.58), (4.59) and (4.60), we conclude that the entropy consistency estimate 
(4.56) holds, provided there exists a positive constant such that 

5 1 N 

(4.61) Ax|| — wjvOM) II 2 + -^\\E f N w N {x,t)\\ 2 D -x < Const. Y A*+$|Auv*i(t)| 2 . 

An energy estimate for the first term on the left is obtained by substituting <Pn{x) = 
j- t w N (x,t) in*(4.45): 

(4.62) ||^» J v(x,()|| ! < 2||^/( U ,„(x,())|| J + 2C„|||- U ,„(x,t)l|5 ) . 

To estimate the second term on the left, we take advantage of its local representation in 
(4.53): since E^Wff(x,t) is a conservative we have 

||^ty^(x,f)|| 2 ,_ 1 = sup || ^^(x,*)^] || 2 

Mz > =1 

( 4<63 ) = sup I (<p{x) - <p N (x), £[f(w N (x, i))]) I 2 < 

ll»»*llo=l 

< (Ax) 2 ||^[/(^(x,i))]||^- l . 

The inequalities (4.61), (4.62) and our assumption (4.57b) imply (4.61), for 

A*IIJfW*r(*.0H 2 + il|^iww(*,OII!)-i < 

< Const. Ax [||/'(u; Ar (x,i)^u; A r(x,i)||5,_ 1 + \\-§^w N {x,t)\\]^ < 

< Const tiL-iv ^ + i|A W „ + i(0| 2 + < 

< Const. '£, 1 2=- N D l , + }.\&w v+ i_{t)\ . □ 

The proof of Theorem 4.6 hinges on estimate (4.61). The latter was verified with the 
help of the first-order estimate 

Q 

\\<p(x) - <p N {x)\\ u < A*||— <p(x)\\ u . 
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In fact, <Pn(x) is a second-order accurate approximation of <p(x), which brings us to the 
final result of this section. 


Theorem 4.7: ( Second-Order Differencing) . The finite-difference/element approxima- 
tion 

(4.64a) iu N (s„_i, t) + ^u N (x u , t) + *)] + 

+ 2AX ( X,/+U 0 ) . — f[ u N{Xv- 1))] = 2 — Qv- > 

is consistent with scalar initial-value problems ( 1.1a ), (1.1b) w.r.t. the quadratic entropy 
U = {| u 2 }, if its viscosity coefficients Q u+ i, satisfy 

1 /•Z„+l 

(4.646) Q v+ 1 > C'on5t.|aJ >+ i (t)| • |Au„ + j(i)| , < + i(i) = [— £ |/"(«jvfc t))| 2 dx] 1/2 . 


Proof: We consider a weak representation of (4.64a) of the form 

(<?*(*)» feM*.*)) + fav(*). &[/(«*(*»*))]) + 

(4.65) 

+ £ N ^ u n{x, t)) Q = 0 , ejv = ^ . 

Here v?(x) = J2k=-N V ^ Xk ^ + ^ Xl+1 ^ Xt+ i (x) I s the piecewise-constant projection of an arbi- 
trary ^-element, ^pat(x) = SitL-jv <£>(£*) A.* (x). 

The local truncation error, En = E^+E^+^-Qn-^, consists of temporal contribution 

(4.66a) ^E$w N (x,t)[tp] = ^?(x) - <p N (x),~w N {x,t)j , 

and spatial contribution E$ = E * N + where 


(4.666) 


£££u>jr(x,t)[<p] = (<p(x) - ^(x), £[/(»*(*.*))]) = 

= (<p{x) -<pW{x),£[f{w N {x,t))} - jzlf{w N {x,t))}) , 


(4.66c) ^^Q N £- v , Jf (x,t)['p] = ^ • 

Once more the entropy dissipation of (4.65) is given by 
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(4.67) 


(§- l w K (x,t),E , i , )w !J (x,t)') + (A w K (z,t),E$wx(x,t)) = 

= |£"- K .D„ + j |A»„ + . (*)P • 

The entropy production does not exceed 
(4.68) 

Jj-ll < w„(x,t), § t E^w N [x,t) > III,-,,,, + II < e H [x,t),E$w H lx,t) >' || t , (l) < 

< <J« 1 ^llfrM*.OI| , + «f||^(*,«)lli,+ i||^ W (x,<)||5,-, , 

and the discrete error is upper bounded by 


(4.69a) 


Now, we note that Q „ + 1 and D u+ ^ = Q„ + i — Q* + i are of the same order of magnitude, 
for by (4.29b), 

(4.696) Q* v+1 _ = /J_ i/2 Q - e 2 ) f"(w u+ i(Z))dt • Au;„ + i(t) < jK + i(t)| • |Auv + i(*)l . 


As before, we choose e n = Const.Ax J. 0 and tjn = Const, with sufficiently large 
Const. > 0, and the entropy consistency (4.56) is reduced to, in view of (4.67), (4.68) 
and (4.69), 


(4.70) A ®|| JU; JV (a;,i)|| 2 + ^||.E&wjr(*»*)||i>-i < Const. £ D v+ i|Aw v+ i(t)| 2 . 

Substituting <Pn{%) = Aio//(x, t) in (4.65), yields an energy-estimate for the first term on 
the left of (4.70); namely, we have 


|| | - t w N {x, t) || 2 < ||| §- t w N (x, t ) || 2 + || £[/( w N (x, t))] - £[f(w N (x, t))J || 2 

+ £ N\l& W N(x,t)l \ 2 + \\§zWn{xMq* 

< \\\$- t W N {x,t)\\ 2 + + Q-ollfewjyr^Olft 7 

which yields, with the help of (4.69b), 


(4.71) || jU^(a;,*)|| 2 < 2{Axf\\^\f{w N {x,t))\\ 2 + \\^w N {x,t)f D 
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The second term on the left of (4.70) does not exceed 
(4-72) 

= sup | (ip{x) -v N {x),^[f{wx(x,t))l - £[/(tu N (i,())l) | ! 
llv>*ll D =1 

Now, since u>n( x , t) is piecewise-linear, we have at each cell, 

a 2 /Aiy v+ i(t)\ 2 

S /"(«*(*,*)) • 

Hence, the inequalities (4.71), (4.72) and our assumption (4.64b) imply (4.70), for 

Asll&wwfo *)!!* + < 

< Const.[ ^-n 1T^{K+ i) 2 I Au W i W I 4 + £ u=-n D u+ i\Aw v+ 1 (t) | 2 ] 

< Const.D^ = _^P l/+ i|Au7„ + i(t)| 2 . □ 

Remarks : 

1. Theorems 4.6 and 4.7 are based on energy-estimating, independently, A ujv(x,t) 
and -§;[f{uN(x,t))). In fact we can do with even less, namely, energy-estimating the sum 
Au^(x,t) + ^[/(ujv(^ji))]' To this end, it will suffice to have a space-time viscosity 
approximation in the “direction” of the local error, A u^(x, t) + Jj[/(ujv(x,t))]. In this 
manner one concludes the consistency of finite-element Streamline-Diffusion (SD) method 
introduced by Hughes and his co-workers and analyzed in [4]. 

2. The case of no mass lumping, (4.27), can be treated similarly. The temporal 

“discretization” in such case is represented by the local truncation error, 

d d 

(4.75a) —E§w n (x, t)[<p] = (tp(x) - (p N { x ), -qI w n( x ’ *)) » 

where <Pn(x) is the ^-interpolant of <p(x) after mass lumping, 

N 

(4.736) (p N (x) = ^2 <p k A. k (x) , (p k = C"V(®ife) • 

k=-N 

The mass lumping in (4.73) adds entropy dissipation to the unlumped mass version in 
(4.64). 

3. Theorem 4.7 verifies the consistency of modern high-resolution nonlinear finite- 
difference approximations, e.g., the TVD and UNO methods described in [2], [3], [14], 
[15], We observe that the viscosity coefficients in (4.64) need not be “limited” at critical 
extrema values. 
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5. INSTABILITY OF THE SCALAR VANISHING VISCOSITY METHOD 


In this section we study the Instability of the scalar vanishing viscosity method. 
To make our point we shall concentrate on the fl^r-framework, where a 7T^-approximate 
solution, ujy(x,t) = Y.k=-N “Jfc(O e,fca: > evolves according to the viscous approximation 


d d d d 

(s.i) + dx^( UN ( x ’ = fo E ” UN ( x ’ ^ + EN dx 


Qn-q^ U n [ x , t) 


In all the examples demonstrated in the previous sections, (5.1) was based on a consistent 
- in fact (at least) first-order accurate approximation of the conservation law (1.1a), 

1 d 

(5.2) || to at (a:) || < Const 0 • . 

In addition, we assume that the viscosity part of (5.1) is parametrized according to the 
setup of Corollaries 4.2 and 4.3. Namely, we have the inequality (4.10) 

—w n (x),Q n —w n (x) \ + ||t0j V (x)|| 2 > Const.£w||— u/Ar(z)|| 2 , 

and the viscosity amplitude, e n, tends “sufficiently slow” to zero 

(5.4) i < ew l 0 . 

In this case, the vanishing viscosity method (5.1) is a consistent approximation of scalar 
initial- value problem (1.1a), (1.1b) w.r.t. all strictly convex entropies, U = Zl a u* This 
assertion was verified in Corollaries 4.2 and 4.3 with the help of the inequality 

(5.5) d N \\-^-w N (x)\\ 2 < f^ tu ^( x )>^ u; ^( :c )) +||wjv(a:)|| 2 , d N = Const.(l - 6)e N . 


Here 6 = < 1 in the general case of Corollary 4.2, see (4.13), and 6 = 0 in the 

entropy conservative case of Corollary 4.3. 

Finally we recall that according to a previous remark, see (3.22), ||iiv(2b £)||l 2 (z) is 
uniformly bounded in time, ||u^(x, < M 2 . Using this together with the fact that 

our ‘preferred’ entropy function is quadratic in the scalar case, U t (u) = |u 2 , enables us to 
upper bound the lefthand-side of (3.15) independently of the mbounds in (3.16), 


pX f Q \ j 

Jt =0 ( g^ U N{ x ,t)iE f N u N {x,t)\ dt = -|Mx,t)||£ a(x) 


t = 0 


l<=r 


< Ml 


Hence, applying (5.5) to twjv(') = ujv(*,t) 


we find after temporal integration that 


(5.6) 


d 1 

- css: 


2 Ml 
(i-o) • 
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The Instability of (5.1) hinges on /Aestimates for the discrete error of the basic 
scheme, 


(5.7) 


1 d 

II II £<*>(*) < Constoo • ^=|| — ttfjv(s)|| , 


and for the discrete error due to the presence of vanishing viscosity, 


(5-8) 


Rn-^w n (x) 


< Constoo • ||tyK(x)||L«>(i) , Rn + Qn = Id-N 


The estimates (5.7) and (5.8) are the strengthened /Aversions of estimates (5.2) and 
(5.3). Indeed, the 2A version of (5.7), 

1 d 

ll^i«itf(*)||i>(.) < Constp- j^\\—w N {x)\\ , 
with p = 2, corresponds to (5.2). The /Aversion of (5.8), 


£ at|| 


dx 


R N j^w n {x) 


lp( x ) ^ Constoo • ||t^iv(^)||Li>(i) , Rn + Qn = Id N , 


with p = 2, yields (5.3) for, 


£n\^w n {x),Qn—Wn{x)] = e:iv||^:^(a:)|| 2 + ^ (^n(s),^: 


dx 


dx 


dx 


Rn-q^w n {x) 


> ejv|| — u>n(z)|| 2 - Constoo • ||wv(x)|| 2 . 


Equipped with these estimates we can iterate on the L p (x) norms of u N (x, t) with the 
help of 

Lemma 5.1: Consider the vanishing viscosity method (5.1) which satisfies the L°°- 
consistency estimates (5.7), (5.8). Then there exists a positive constant (independent of N 
and p), Const. > 0, such that for any even integer p > 2 we have 


(5.9) 


Vp 


I|«nOM)IU*(*) < e ConstA |Hu w (*,t = 0)||x, w + Const. -j==—£ 

Proof: Multiplying (5.1) by pu p I f 1 (x,t) and integrating by parts we obtain 

4lMMIIL(.) = -p[p-l)S2=oU P N 2 {x,t)§;UN(x,t)E f N u N {x,t)dx 


(5.10) 


-p(p - l)e w fH o u p N 2 (x, t) \j^u N (x, t)f dx 
P^N /i=o ^ _1 (x, t)£; [R N §^u N (x,t)] dx 


= I + 11+ III. 
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For even integers p> 2, the weighted Cauchy-Schwartz inequality is used to upper bound 
the first term on the right 


(5.11) 

I < Pip- 1 ) £ n J x _ o u P n 2 (x, t ) 


d t . 

— u N {x,t) 


dx + E^E — 11 f U P N 2 ( Zj t ) • [E f N u N {x, t)fdx , 
Jx= 0 


and together with the second term on the right we have, in view of (5.7), 

(5.12) 1 + II < Constl • Wlv-Hx) ' II *)!&(«) * 

Using (5.8) we can Holder the third term, 


(5.13) 


III < pe N \\u p N 1 (ar,*)ll Xr5 * r(x) llfe [**&«*(*,*)] llx»(x) ^ 

< Cons^ • p • ||«^(s, <)||£ P(l) . 

Inserting (5.12) and (5.13) into (5.10) we obtain, after division by a common factor of 

SIM*.*)II5?.), 

d v d 

( 5 - 14 ) m\Wn[xMIv{z) < Cousin • -j^\\—u N {x>t)\$fi {x) + 2Const 00 \\u N {x,t)\\ 2 LT[x) . 
Temporal integration of the last inequality yields 


2Constoo-t 


• f ||utf(x,t)||L- (z) < II u N {x,t = 0 )||£ p(l) + 

+Const\ Q • e~ 2Con,tc °' T £ N • ||^u^(z,r)||* 2(;i;) dr 

and (5.9) follows with the help of (5.6). □ 


Lemma 5.1 shows that the L p (z)-norms of u^(x,t) are bounded (w.r.t. N,p and t ), 
at least for “sufficiently high” L p (a:) -norms, provided ejy tends to zero “sufficiently slow”. 
Specifically, if instead of (5.4) we have the stronger 

(5.15) ^ e N i 0 * 
then ||« 7 v(x,t)||£i.(j : ) are bounded for p < log N, 

(5.16) ||ujv(x,f)||ip( x ) < M(t ) , p < log IV . 

Now, we assert that the L°°(x)-norm of 7TAr-polynomials does not exceed a constant, say 
10, times their “sufficiently high” L p (x)-norm, say p ~ log N, for by one of the Gagliardo- 
Nirenberg inequalities, e.g., [16, Section 3] 
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< («?)* ■ IIAM*)II% • IKMIlgf.) < 

< l.l-JV*||«, N (x)||* ) -|K( a; )||£ff, ) < 

< 1.1 • e 2 ||wjv( x )IU»’= , o* w (x) . 

We conclude 

Theorem 5.2: (L°°- stability) .// the scalar vanishing viscosity method (5.1) satisfies 
(5.7), (5.8) and (5.15), then it is L°°-stable, i.e., there exists a constant (independent of 
N), Mm, such that 

||ti*(z, OIU^M < Moo . 


Remarks : 

1. The consistency estimate (5.7) together with (5.6) imply 

(5.18) f \\Ej f u N (x i t)\\la B i x) dt < Const. -J— . 

Jt = 0 v ' 1 V6n 

Thus, the presence of vanishing viscosity in (5.4), parametrized according to (5.15), guar- 
antees that the local error of the basic approximation in (5.1) tends to zero. 


j ^_ q \\En u n(x, t) \\L°°[ x )dt < Const. 


0 . 


y/m^N 

2. If instead of (5.15), the viscosity amplitude e n is restricted by the weaker (5.4), then 
(5.18) yields 


J t _ ||-E&ttw(s,t ) Wlootfdt < Const. , 

and we conjecture that the vanishing viscosity method (5.1) remains L°°-stable in such 
case. 


The first step in implementing Theorem 5.2 requires us to verify the L°°-consistency 
estimate (5.7). To this end let us consider basic discretization operators which operate 
linearly on the flux /(•), 

(5.19) P^w N {x) = Pn f (wn (®) ) 

We recall that the spectral, ^dospectral, finite-difference and element methods discussed 
in Section 1, are the canonical examples of basic approximations which belong to this 
“linear” category. The consistency estimate(5.7) then reads 

1 d 

\\[I - P N }f{w N (x))\\ L ^( x ) < Constoo • ^=!I^M*)|| , 
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with a constant Const ^ which may depend on the L°°-bound on w N (x) but otherwise is 
dependent on N. Consequently, estimate (5.7) is fulfilled provided for all 10 ( 2 ) in JT 1 [0, 27 r] 
we have 

(5.20) ||[/-P^]u;( 2 )||i,oo( s ) < Const. -~\\-^-w(x)\\ . 

We shall verify the consistency estimate (5.7) - or equivalently (5.20), using the following 
von Newmann like stability analysis. To this end we define the symbols 

(5.21) g Nik {x) = e~ ikx P N e ikx , 
and state 

Lemma 5.3: The linear basic approximation = P^/(*), fulfills the L°° -consistency 
estimate (5.7), provided its symbols satisfy 

(5.22) max|l — <7w,jfe(s)| < Const, min ^1, , |/e| < 00 . 

Proof: For |&| > N we have by (5.22), 

II £|*i>* M k ) • [! - 9NA x )\ eikx \\\°°{x) < 

(5.23a) < 2T|jfc|>7v A; 2 |u/(A;) | 2 • |1 — gN,k{x)\ 2 • £|*|>w pr < 

< *^£|*|>** 2 NA0| 2 . 

Moreover, (5.22) implies that for |&| < N we have 

II E w <„ *(*) • [1 - < N • E w <„ |A(*)|* • |i - < 

(5.236) 

< S 2 f ti E W < W * ! |&(*)| 2 • 

Hence (5.20) follows in view of (5.23a), (5.23b), 

||(J - Pit) w(*)Hio. w = II ££L-oo H k ) • [1 - ^(*)]«“*lli» w < 

< **K*)I* = [2Const.^\\fw{x)\\\ 2 . 

□ 

We observe that the exponentials {e ,A:z }iL-Ar are eigenfunctions of the linear discretiza- 
tions operators associated with the spectral method, P^r = S jy, V’dospectral method, 
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P N = %l) N , finite-difference method, P h = FD^ X and finite-element method (based on 
quadratic “preferred” entropy t/*(u) = |u 2 ), Pat = FE& X . We claim that therefore (5.22) 
is satisfied, which in turn implies that the L°°-consistency estimate (5.7) is valid in these 
cases. Indeed, since P ^ are bounded 


(5.24a) max|< 7 jy,k(a:)| = max \Pse ,kx \ < Const. 

Moreover, by the first-order accuracy estimate (5.2) we have 

Ill'll - »w,»(*)llli>(.) = l |[/- ■ 

Now, if are eigenfunctions of Pat, then gu,k(^) are constant amplification factors, 

9N,k{z) = 9N,ky and the last inequality yields 

1 k 2 

(5.246) max |1 - gN.k{ x )\ 2 < —Constl • , \k\ < N . 

This together with (5.24a) imply (5.22) as asserted. 


Example 5.4: L^-stability of the scalar SV method . We consider the SV method 
(4.17), (4.23) 

d d \ d d 

(5.25a) — [t*jv(a;, i)l + ^[JV/(«w(®,*))l = [Qjv(a=,t) * ^u N {x,t)} , 


where Qn(x, t) = 2|k|= mN Qk[t) e ' kx is a viscosity kernel with monotonically increasing 
coefficients such that 


(5.256) Qk[t) > Const. - - -r —. — — , fc 2 > m 2 N = — — 

v k 2 • log JV “ N Const, log N 




The basic spectral or i/>dospectral approximation satisfies the L°°-consistency estimate 
(5.7) by Lemma 5.3. Next, we verify estimate (5.8) with Rn + Qs = Const.Id n. In this 
case, the corresponding kernel R^(x,t) is given by 


Rn{x , t) = Const. e ' kx + X) Rk(t)e ,kx = R N {x) + R^{x, t) , 

\k\<m N |*|>m N 

where Pjy(x) is a multiple of the 7r mjv -Dirichlet kernel, 


R n (x) = Const. D mN (s) , 

and by (5.25b), R^(x,t) has monotonically decreasing coefficients 

- NP 

° - 1 Rk W - k 2 -lo g N ’ ^ ~ mN ■ 
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Since iE^x) * wn{%) is a 7T mw -polynomial, we can estimate its derivatives by [23] 

e\\h\ R N{. x ) * ^ e N m N ’ Pjv(*)IU*(*) • ||w JV (x)||l"o( i) < 

(5.26a) < Const.e^m^logm^Ijiu^^JIjioo^) < 

< Constoo ||u;jv(x) |U°°(*) • 

Since the coefficients of logiV • N~^f^R'}f(x,t) are monotonically decreasing, 

0 jlog N‘N~ fi k 2 R k (t) < 1 , 
we can apply [10, Lemma A.l] to obtain 

MsiWfoO * fe«»Jr(*)]|li».(*) < 

< ^££77 • li log AT • N-P£rR+{x,t)\\ LHx) • ||u;^(z)|| i « (l) < 

(5.266) 

< • Const, log N • ||tajv(x)|| L ~ (:c ) < 

< Const oo • ||taAr(*)|U~>(it) • 

Thus, estimate (5.8) holds in view of (5.26a), (5.26b). Finally, since (3 < the viscosity 

amplitude — AT -/3 satisfies (5.15), and Theorem 5.2 guarantees the L°°-stability of the 
SV method (5.25a), (5.25b). 
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6. CONSISTENCY AND L°°-STABILITY IMPLY CONVERGENCE 
We begin with 

Theorem 6.1: Consider the semi-discrete approximation (1.4). We assume that 

(i) The approximation (1.4) is consistent with the initial-value problem (1.1a), (1.1b) w.r.t. a 
“ relevant ” class of entropy functions, U, 

and 

(ii) The approximation (1.4) is L°° -stable. 

Then for any C 2 entropy pair ( U,F ) associated with (l.l), such that U (E £/, 

d d 

(6.1) — [t/(u N (x,t)j + — [F(utf(x,t))] , U eU , 

belongs to a compact subset of Hf 0 \[x,t). 

Proof: Multiplying (1.4) by t/'(u//(x,f)) we obtain (3.1), which we rewrite us 
&[If(u w (x,i))] + £[F(ujy(*,i))] = 

( 6 - 2 ) = h < U'{u N (x,t)),E f N u N (x,t) > - < \u N {x,t)),E f N u N (x,t) >= 

= I + 11. 

By the first part of Lemma 3.3 we have, see (3.9), 


II I — dx < U'{ u N[ x it))i J^N u N(x,t) > ||tf-i(£jt) < 

# . IOC 

(6.3) 

< ||17'(ujv(a:, t))||z,“ ( x ,t) * Consti • y/efj — > 0 

and hence the term I lies in the compact Hf 0 l. The third part of Lemma 3.3 gives us 


(6.4) || II = - < — [U'(u N {x,t))},E f N u N [x,t) > \\ L } oc {x,t) < Constni , 

and this estimate implies, with the help of Murat’s Lemma [22], that the term II also lies 
in the compact of Hf 0 ], which completes the proof. □ 

One can use now compensated compactness arguments [22], in order to turn the con- 
clusion of Theorem 6.1 into a convergence proof. The current state of the art of these 
arguments in this context, includes scalar and 2x2 systems of conservation laws [22], [l]. 
We start with the scalar case. 


Theorem 6.2: Consider an L°°-stable approximation (1.4) which is consistent with the 
scalar initial-value problem (1.1a), (1.1b) w.r.t. all strictly convex entropies, U = U a \\. Then 
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(a subsequence of) u^(x,t) converges weakly to a weak solution, u(x,t), of the conservation 
law (l.la), 

(6.5) ^[«(*, 0)1 = °- 

Proof: Since the right-hand sides of (1.4), (3.1) were shown to lie in the compact of 
Hf 0 \ (x, t ), we can apply the div-and lemma [22] to their left-hand sides: abbreviating g(u) = 
u)Yim g{uff), then for any C 2 entropy pair ( U(u),F(u )) we can extract a subsequence (still 
denoted by u^) such that 

(6.6) uF(u ) — U(u)f(u) = u • F{u) — U(u) ■ f(u) 

But (6.6) depends continuously on ( U , F) in the C 1 topology and therefore it remains valid 
for piecewise C x entropy pairs. Following Krushkov [6] we choose 

U{u ) = |u — c| , F(u) = sgn (u — c) • ( f(u ) — f{u)) c = Const , 

in which case (6.6) reads 

u • sgn (u - c) • (/(u) - /(c)) - |u — c| • /(u) 

= u • sgn (u - c) • ( f(u ) - /(c)) - |u - c| • f{u) . 

Equivalently, we can rewrite this as 


(u - u)sgn (u - c) ■ (/(u) - /(c)) = |u - c| • (/(u) - f(u )) . 

Let us examine the last equality by restricting the weak limits to an arbitrary fixed (x, t) 
location; with c = u(x, t ) we find after little rearrangement that 

|u - u| • ( f(u ) - f(u)) = 0 . 

This implies that 

(6.7) f{u) = /(u) , 

for otherwise |u — u|(x, t ) = 0, which in turn leads again to (6.7) 3 Taking the weak limit of 
(1.4), then (3.9) implies that the righthand-side tends to zero, and by (6.7) the lefthand-side 
amounts to having (6.5) , 

j^[u(x, t)] + ~[f(u{x, *))] = j^[«(*, t)] + ^~[f(u(x, i))] = 0 . □ 

3 If |u — u|(i, t) = 0 then y(u)(x, t) = g(u(x, t)) for any C 1 function. This follows, for example, by noting 
that the associated probability measure is concentrated at the single point u, consult [22]. 
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The above compensated compactness argument for convergence is due to Tartar [22, 
Theorem 2.6] who made use of (6.6) to further deduce that in case the scalar flux /(«), is 
nonlinear, i.e., when there is no interval on which f'(u) = Const., then the convergence 
Un{x, t) — * v(x,t) is in fact strong in L* oc (x,t), p < oo. Moreover, if the scalar flux /(«) 
is genuinely nonlinear (GNL), i.e., f"(u) ^ 0, then the same conclusion of strong L\ oc [x,t) 
convergence holds, by applying the div-curl lemma to a single strictly convex entropy 
function U(u) in U [22, Remark 30]. 

Theorem 6.3: Consider an L°°-stable approximation (1.4) which is consistent with the 
nonlinear scalar initial value problem (1.1a), (1.1b) w.r.t. a nonempty u relevant ” class of 
strictly convex entropy functions U. Then we have 

A. Convergence : If either U contains all strictly convex entropies, U = Zi a ii or f(u) is GNL, 
then (a subsequence of) ujv(x, t) converges strongly in Lf oc (x,t), p < oo, to a weak solution 
u(x,t) of the conservation law (1.1a), 

(6.8) y t [u{x , *)] + j^[/(u(x , t ))] = 0 . 

B. Entropy inequality : For each entropy function, U[u ) in U, whose entropy production 
tends weakly to zero 

' Q 

(6.9 a) w lim - < — U'{u N {x,t)) , E f N u N [x,t) >~ = 0 , 

iV— ♦ OO ox 

the weak limit solution u[x, t) satisfies the entropy inequality 

(6.96) S iC7( “ (l ’ i))] + !)) 1 - ° ' 

Proof: The amount of entropy dissipated by (1.4), 

d 

- < Q^U'(u N (x,t)) , E f N u N (x,t) > + 

tends weakly to a negative measure in view of the L} oe (x,t) bound in (3.11). Adding this 
to (6.9a) we conclude that the second term on the right of (6.2) , tends weakly to a negative 
measure, 

Q 

(6.10) wlim^ln = - < —U'{u N (x, t)) , E f N u N (x,t) >] < 0 . 

Thus, in view of (6.3), (6.10) and the strong convergence, u N (x,t ) — > u(x,t), the weak 
limit of (6.2) recovers the entropy inequality (6.9b). □ 
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Example 6.4: Convergence of the SV methods. We consider the SV method (4.17), 
(4.23) 

d d \ d d 1 

(6. Ha) qi\un{x, <)] + ~[P N f(u N (x, *)] = -q^{Qn{x, t ) * q^un(x, *)] , fi < - , 

where Qn(x, t) = X)|jt|= mjv Qk[t)e tkx is a viscosity kernel with monotonically increasing 
coefficients, 

(6.116) 

- N p N p 

Qk{t) > Const. - -r-r— — , k 2 >m 2 N = — — — , Const. > Const 0 . 

k* • log N Const, log N 

The SV method (6.11a), (6.11b) satisfies 

(i) It is consistent with the scalar initial-value (l.la), (l.lb) w.r.t. all strictly convex en- 
tropies (by Theorem 4.3), 

(ii) It is L°°-stable (by Theorem 5.2), 

and as a consequence of (i) ,(ii) and the first part of Theorem 6.3, 

(iii) The SV solution, un(x, t), converges strongly to a weak solution, u(x,t), of the non- 
linear scalar conservation law (l.la). 

(iv) The entropy produced by the basic approximation of (6.11a) tends strongly to zero 
(by Theorem 4.3 III). 

(v) The entropy produced by the viscosity approximation of (6.11a) tends weakly to zero, 
consult [10, Section 5]. 

As a consequence of (iii) , (iv) , (v) and the second part of Theorem 6.3 we conclude: 
The SV method (6.11a), (6.116) converges strongly to the unique entropy solution of the 
nonlinear scalar conservation law (l.l). 

Remarks : 

1. The last conclusion extends the convergence results of [20], [10] which was restricted 
to the inviscid Burgers’ equation where f(u) = |u 2 . 

2. The instability of the SV method (6.11a), (6.11b) with (3 < 1, would imply its 
convergence along the lines of Example 6.4. 

Example 6.5: Convergence of Centered Finite-Difference Methods . We consider the 
conservative finite-difference/element method (4.27) 

■^U^{x i/,t) + 2Ax'[/( li iv( a 'i/+l> 0) /(^n(Xi/-1) t))] = 

(6.12a) 

= 2A*[ < 2«'+j^ U *'+jW ~Qv -^ u J 
with numerical viscosity coefficients, Q v+ l < Qooi such that 
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( 6 . 126 ) 


<?„ +i > qu i + Cmst - . <32+1 = 2 e/*(«wj(e)W • 

This method is consistent with the scalar initial-value problem (1.1a, (1.1b) w.r.t. all strictly 
convex entropies (by Theorem 4.5). Since (6.12a), (6.12b) is also known to be L°°-stable, 
e.g., [14], the first part of Theorem 6.3 implies strong Lf oc (x, t) convergence to a weak limit 
solution, u(x,t) of (1.1a). Finally, since (6.11a), (6.11b) contains more numerical viscosity 
than the entropy conservative schemes (1.24b), then u(x,t) satisfies (6.9b), (consult [18] 
21]). We conclude: 

The finite-difference method (6.12a), (6.12b) converges strongly to the unique entropy 
solution of the nonlinear scalar conservation law (l.l). 

Example 6.6: Convergence of Upwind Differencing . We consider the conservative 
finite-difference/element method 

[g u N( a 'X'-lj 0 H" gUjy(x !/,t) + gU;y(x I/ -(.i, t)j + 

(6.13a) 

"^2Azt/( U 3v(^+lj t)) — /(ujv(x v _x, t)] = 2 Az |^ U v+i(^) — Qv- iAu„_i(t)] , 
with numerical viscosity coefficients, Q u+ 1 < Q^, such that 

r l f x *+i „ 1^/2 

(6.136) Q„ + i > C'onsf.|a I , + i| , a„ + i = ^ \f\u N {x,t))\ 2 dx 

This method is consistent with the scalar initial-value problem (l.la), (l.lb) w.r.t. the 
quadratic entropy U = {|u 2 } (by Theorem 4.6). Since (6.13a), (6.13b) is known to be In- 
stable, e.g., [14], the first part of Theorem 6.3 implies strong .Lf oc (x, t ) convergence to a weak 
limit solution, u(x, t), of the nonlinear conservation law (l.la). Moreover, the quadratic 
entropy produced by the basic approximation of (6.13a) tends strongly (in H^l(x, t)) to 
zero, while the viscosity approximation is purely dissipative, consult Theorem 4.6. Hence 
with the help of the second part of Theorem 6.3 we conclude: 

The finite-difference/element method (6.13a), (6.13b) converges strongly to the unique 
entropy solution of the GNL scalar conservation law (l.l). 

Example 6.7: Convergence of Second-Order Differencing. We consider the conserva- 
tive finite-difference/ element method 


dt [e 6 Wjv(X|/,t) + gU^Xy+i, t)j + 

+ ih[f^N{x v+ 1 ,t)) - /(Utf(x„_i,<)] = 2^[Q„ + iAu„ + i(t) - Q^iAu^dt)] , 


51 




with numerical viscosity coefficients, Q v+ l < <5oo> such that 


(6.146) Q v+ i > Const. \a? v+ i_\ • |Aw„ + i (i)| , 


a 







This method is consistent with the scalar initial-value problem (1.1a), (l.lb) w.r.t. the 
quadratic entropy U = {|u 2 } (by Theorem 4.7). Now, if (6.14a), (6.14b) is L°°-stable, 
then the first part of Theorem 6.3 implies strong Lf oc (x, t) convergence to a weak limit 
solution, u(x,t), of the nonlinear conservation law (l.la). Moreover, the quadratic entropy 
produced by the basic approximation of (6.14a) tends strongly (in Hf 0 \{x, t)) to zero, while 
the viscosity approximation is purely dissipative, consult Theorem 4.7. Hence with the help 
of the second part of Theorem 6.3 we conclude: 

The finite-difference/element method (6.14a), (6.14b) converges strongly to the unique 
entropy solution of the GNL scalar conservation law (1.1), provided it is L°°-stable. 


Remarks : 

1. The conclusions of Examples 6.6 and 6.7 remain valid with or without mass lumping 
on the left of (6.13a), (6.14a). 

2. The question of L°°-stability for the second-order accurate approximation (6.14) 
remains open. 


We conclude this section with convergence results for strictly hyperbolic GNL systems 
of two conservation laws. Making use of Theorem 6.1 and DiPerna’s results in [1] we have 

Theorem 6.8: Consider an L°°-stable approximation (1.4) which is consistent with 
the 2x2 initial-value problem (l.la), (1.16), w.r.t. a nonempty “relevant” class of strictly 
convex entropy function U. Then we have 

I. If ||t£w(x, i) ||x“» c (iE.i) is sufficiently small then u^{x,t) converges a.e. to an admissible 
(i.e., entropy satisfying) solution, u(x, t), of (l.la). 

II. If (l.la) is equipped with quasi-convex Riemann invariants and (5.6) holds, then 
Ujy(x, t) converges a.e. to an admissible solution u(x,t), of {l.la). 

Examples of utilizing Theorem 6.8 in the context of the vanishing viscosity method for 
2x2 systems will be given elsewhere; the consistency analysis in such cases can be carried 
out more conveniently in terms of the entropy variables, v, instead of the conservative 
variables, u [21 j. 
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7. APPENDIX: CONSISTENCY WITH THE CONSERVATION LAW 


In this section we study the order of accuracy for various types of discretizations, in- 
cluding spectral and pseudospectral methods, as well as arbitrary recipes of finite-difference 
and finite-element methods which admit the general conservative form (1.14), (1.15). It 
is shown that all these methods are at least first-order accurate and hence consistent 
approximations of the conservation law (1.1a). 

We start with 

Theorem 7.1: The spectral method (1.10), (1.11) is infinitely-order accurate, i.e., 
accurate of any order s > 0. 

Proof: With fk(t) denoting the k-th Fourier coefficient of /(ty^(x)), we have for any 
integer r > 0, 

r>_r oo 1 3r 

(7.1) ||(J-s„)/(«. K M)|| J = 2,r ^ |A| ! <^ £ l*f r IAI ! = 4;ll£;/(«»W)ll ! . 

\k\>N iV k— — oo iV OX 

and (1.6) follows with C 0 = ||/ , (^(®))||l“(z}- □ 

Remark: The linear case f(u) = u shows that the above accuracy estimate is indeed 
sharp. We note that this estimate makes use of the L°°-bound of u^(x, t ). In the quadratic 
scalar case, f(u) = |u 2 , a weaker accuracy estimate of order s = 1 was proved in [20, 
Le mm a 3.1], independently of L°°-bounds. 

We turn to the pseudospectral method (1.12), (1.13) where P^vjn(x) = 0w/(ti;Ar(x)) 
is the ^-interpolation of /(u>;v(x)) at x u+ $ = [y + 0)Ax, u — 0,1, •••,2 N with fixed 
0 < 9 < 1. We recall that 

(7.2a) v = SV + -A at, 

where the aliasing projection, A #, is given by 

(7.26) A n w{x) = E (E *(* + j(2N + l))e ii2 * 9 ]e ikx . 

|i|<AT jyo 

To treat the pseudospectral case, we first prepare 

Lemma 7.2: There exist constants (independent of N), C r , such that for any r > | 
we have 

(7-3) \\A N w(x)\\<C r ~\\(I- S N )-^w[x)\\, r>i 
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Proof: Using the aliasing relation (7.2b) we have 


II4*®(*)IP = E W <* I + J( w + l)e« J 'T < 

(7.4) 

< E WS »E*» • E//o |/P'JV !r |t»(* + j(2N + 1)P], 

Since \j\N < \k + j(2N + l)\ for \k\ < N,j 7 = 0, the right-hand side of (7.4) does not exceed 
E 171375757 £ El*+j( 2 ^+l)P r N(* + i( 2 ^+l)| ! = 

jyo mi iV |*|<jvjyo 

= ^-EtL- e pp>mp. 

iv jyo I-/ 1 |«|>iV 

and (7.3) follows with C r = (X^o 

Equipped with Lemma 7.2, we can show 

Theorem 7.3: The pseudospectral method (1.12), (IAS) is “infinitely-order” accurate, 
i.e., accurate of any order s > |. 

Proof: Since and I — Sn are orthogonal projections, we have 
(7- 5 ) ||(/-^)/(tw jV (a:))|| a = ||A^/(u; iv (x))|| 2 -|- ||(J- 5^)/(u;^(x))|| 2 ; 

by Lemma 7.2, the first term on the right does not exceed 

||A„/(«.„(x))|p < Cl ■ j£||(7-S W )£r/(»«W)|p < 

C r^ll^/(^W)lP. 

and together with upper bound of the second term from Theorem 7.1, the accuracy estimate 
(1.6) follows with C 0 = (1 + C' r )||/ , (^(x))j| i co (l ).n 

Next, we consider general conservative discretizations of the form (1.14), (1.15). These 
approximations can be interpreted as evolution schemes in n N governed by (1.3), with 
which is defined according to (1.17), 


(7.6a) 


N ■, 2N 

pfa»(z) = £ *»«“*. h = E\K + i - K-d‘~ a *‘- 


k =-N 


2nik " + 5 ' 2 


Here h f v+ , = h f {wjf{x v - p+ \), • • • ,w N {x v+p )) is the Lipschitz continuous numerical flux 


which is consistent with the differential one 

(7.66) h f {w N {x v ),---,w N (xJ)) = f{w N {x„)). 
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We claim that such schemes are at least first-order accurate with (1.1a). To see this we 
first prepare 


Lemma 7.4: Let 
(7.7) 


i"M = “(fMf. A* 


2tt 


2N + l’ 

denotes the sliding average of a 2n-periodic function tv(x). Then we have 
(7.8a) ||tF(x)|| < ||u;(x)||, 


(7.86) 


M*) - u/(x)|| < Const. E|| |^||, 


r < 2. 


Proof: Using the obvious inequality (7.8a) we have 

(7.9) ||u>(x) _ w(x)|| ^ 2||io(x) _ <SWio(x)|| + ||S' w io(x) - 3Sf5T(*)||. 

By Theorem 7.1, the first term on the right is bounded by 777-II §^||* Also, for w^(x) = 
Snw(x) we have 

(7.10) w N (x) = E &(/;) 

o<1*|<jv 


sinker 


N 


2 . Akx 


fc- Al 


w N {x) = E ™{k)e tks -, 


2 k~-N 

hence the second term on the right of (7.9) does not exceed 

cir) 1 

(7.11) |K(*)-o»(x)||< E l*(*)l , |i- = iie-l , <5(A*) 4 i:w>(*)P. 

0<|jfc|<JV K 2 D k 

and (7.8b) follows with Const = 2+^-.D 
We arrive at 

Theorem 7.5: The conservative approximation (1.14), (1-15) ,s least first-order 
accurate with the conservation law (1.1a). 

Proof: Let us recall (1.16a) where JEf^(x) denotes the JV-trigonometric polynomial 
which interpolates H^{x„) = h f 1, i.e., 


(7.12) 

We then have 

(7.13) 


N 

£ 

k=—N 


A n- 2iV 

H N (x) = E H k e*\ H k = ^ E ^ 

Z7r *=0 2 


H/OMx)) - P^(x)|| < J + 7/ + III + IV. 
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By Theorem 7.1 we have 


(7.14) I = ||/(tMx)) ~ Sjv/(u> w (*))ll < ^11 , r < 2 . 

Lemma 7.4 implies 

(7.15) 11= ||Sjv/(ttfjv(x)) - S N f(w N (x))\\ < Const.-^\\-^f(w N (x))\\ , r < 2 . 

Let rp N f(w N (x)) denotes the 7r^-interpolant of /(u/jy(x)) at x„ + i = (u + |)Ax, v = 

0, 1, • • • , 2 N. Lemma 7.2 and 7.4 imply 

(7-16) 

III = |6V/(wn(x))-^ iV /(ui a -(x))|| < ||A.v/fw iV (x))|[ < Const.-^\\~f(w fJ (x))\\ , r < 2 

Finally, we recall (1.16b) asserting that Pj^w N {x) is nothing but the sliding average of 
H n (x), i.e., 

(7.17) PffW N {x) = Hn{x). 

Consequently, for the fourth term on the right of (7.13) we have 

IV = || rp N f(w N (x)) -PffW N (x)\\ < \\tp N f{w N (x)) -J^(*)|| = 

( 7 * 18 ) 

= [£ 2 =o \f{w N {x v+ i)) - h f v+1 _ I 2 Ax]*. 

Using the Lipshitz continuity of the consistent flux h 1 i we can upper bound 
\f{ w N(z v+ i))-h f v+L \ 2 < L 2 [J2 \w N (x v+t )~w N (x v )\] 2 <2pL 2 - \wn{x v -h)-vn{x v )\ 2 , 

3 t=-p t=- P 

and hence the right-hand side of (7.18) does not exceed 

p-i 2 N q 

(7.19) IV < (2p ■ L 2 • J2 |u>jv(zj/+<) “ «Mxp)| 2 Ax)s < 2 pL • Const. Ax||— iujv(x)||, 

0=—p u—0 ° X 

and the first-order accuracy estimate follows from (7.13) in view of (7.14), (7.15), (7.16), 
and (7.19). □ 

We observe that the first three terms on the right of (7.13) are in fact second order accu- 
rate, and hence the second-order accuracy of the conservative approximation (1.14), (1.15) 
depends on the ^-distance, (7.18), between the numerical flus h^ +L , and the midvalues 
f(w N (x l/+ i)). Making use of the detailed structure of in (4.26) we conclude 
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Corollary 7.6: The conservative approximation (4.27) is second-order accurate with 
the conservation law (1.1a), provided its viscosity coefficients satisfy 

(7.20) ||Q v+ i|| < Const. | Atw^il , Q„ + i = Q u+ i{w n {x),x) . 

Proof: Inserting (4.26) into (7.18) we obtain 
(7.21a) 

IV = [E^ol/MaW*)) “ \ ifi.w N (x u )) + f{w N [x u+1 ))) + \Q„+ iAw w+ i_\ 2 Ax} 1/2 < 

< (Ax) 2 [Const.ll^H + (E^= 0 |^tQ v+ ^Au; v+ i| 2 Ax) 1/2 ] . 

Assumption (7.20) implies that the righthand-side of (7.21a) does not exceed 

(7-216) C 'o*||^2« , Jv(a:)|| + C , 1 .||^ti; JV (x)|| , 

and the result follows from (7.13) in view of (7.14), (7.15), (7.16) and (7.21). □ 

The essentially nonlinear high-resolution approximations surveyed in, e.g., [3], [15], are 
characterized by viscosity coefficients, Q u+ 1 , which satisfy (7.21) at all but a finite number 
of gridcells. In this case we have, in agreement with [3], 

Corollary 7.7: The high resolution approximations (4.27), whose viscosity coefficients 
Q v+ 1 satisfy (7.21) at all but a finite number of critical gridcells, are accurate of order 
s = 3/2 with the conservation law (1.1a). 

Verification of this corollary is immediate noting that the contribution of the finite 
number of critical gridcells to the summation in (7.21a) is of order (Ax) 3 / 2 . 
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